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Abstract 

We review the formalism of the effective average action in quantum field theory 
which corresponds to a coarse grained free energy in statistical mechanics. The as- 
sociated exact renormalization group equation and possible nonperturbative approx- 
imations for its solution are discussed. This is applied to QCD where one observes 
the consecutive emergence of mesonic bound states and spontaneous chiral symmetry 
breaking as the coarse graining scale is lowered. We finally present a study of the 
chiral phase transition in two flavor QCD. A precision estimate of the universal crit- 
ical equation of state for the three-dimensional 0(4) Heisenberg model is presented. 
We explicitly connect the 0(4) universal behavior near the critical temperature and 
zero quark mass with the physics at zero temperature and a realistic pion mass. For 
realistic quark masses the pion correlation length near T c turns out to be smaller 
than its zero temperature value. 
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1 Effective average action 



Quantum chromodynamics (QCD) describes qualitatively different physics at different 
length scales. At short distances the relevant degrees of freedom are quarks and glu- 
ons which can be treated perturbatively. At long distances we observe hadrons, and an 
essential part of the dynamics can be encoded in the masses and interactions of mesons. 
Any attempt to deal with this situation analytically and to predict the meson properties 
from the short distance physics (as functions of the strong gauge coupling a s and the cur- 
rent quark masses m q ) has to bridge the gap between two qualitatively different effective 
descriptions. Two basic problems have to be mastered for an extrapolation from short 
distance QCD to mesonic length scales: 

• The effective couplings change with scale. This does not only concern the running 
gauge coupling, but also the coefficients of non-renormalizable operators as, for ex- 
ample, four quark operators. Typically, these non-renormalizable terms become 
important in the momentum range where a s is strong and deviate substantially from 
their perturbative values. Consider the four-point function which results after inte- 
grating out the gluons. For heavy quarks it contains the information about the shape 
of the heavy quark potential whereas for light quarks the complicated spectrum of 
light mesons and chiral symmetry breaking are encoded in it. At distance scales 
around lfm one expects that the effective action resembles very little the form of the 
classical QCD action which is relevant at short distances. 

• Not only the couplings, but even the relevant variables or degrees of freedom are 
different for long distance and short distance QCD. It seems forbiddingly difficult 
to describe the low-energy scattering of two mesons in a language of quarks and 
gluons only. An appropriate analytical field theoretical method should be capable of 
introducing field variables for composite objects such as mesons. 

A conceptually very appealing idea for our task is the block-spin action ^ It realizes 
that physics with a given characteristic length scale I is conveniently described by a func- 
tional integral with an ultraviolet (UV) cutoff A for the momenta. Here A should be larger 
than l~ l but not necessarily by a large factor. The Wilsonian effective action replaces 
then the classical action in the functional integral. It is obtained by integrating out the 
fluctuations with momenta q 2 ^ A 2 . An exact renormalization group equation 0, [|, |], [5], || 
describes how changes with the UV cutoff A. 

We will use here the somewhat different but related concept of the effective average 
action @ which, in the language of statistical physics, is a coarse grained free energy 
with coarse graining scale k. The effective average action is based on the quantum field 
theoretical concept of the effective action || T which is obtained by integrating out all 
quantum fluctuations. The effective action contains all information about masses, cou- 
plings, form factors and so on, since it is the generating functional of the 1PI Green 
functions. The field equations derived from T are exact including all quantum effects. For 
a field theoretical description of thermal equilibrium this concept is easily generalized to a 
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temperature dependent effective action which includes now also the thermal fluctuations. 
In statistical physics T describes the free energy as a functional of some convenient (space 
dependent) order parameter, for instance the magnetization. In particular, the behavior 
of T for a constant order parameter (the effective potential) specifies the equation of state. 
The effective average action Tk is a simple generalization of the effective action, with the 
distinction that only quantum fluctuations with momenta q 2 <: k 2 are included. This 
can be achieved by introducing an explicit infrared cutoff ~ k in the functional integral 
defining the partition function (or the generating functional for the n-point functions). 
Typically, this IR-cutoff is quadratic in the fields and modifies the inverse propagator, 
for example by adding a mass-like term ~ k 2 . The effective average action can then be 
defined in complete analogy to the effective action (via a Legendre transformation of the 
logarithm of the partition function). The mass-like term in the propagator suppresses the 
contributions from the small momentum modes with q 2 ^ k 2 and Tj, accounts effectively 
only for the fluctuations with q 2 ^ k 2 . 

Following the behavior of T k for different k is like looking at the world through a 
microscope with variable resolution: For large k one has a very precise resolution ~ k^ 1 
but one also studies effectively only a small volume ~ k d . Taking in QCD the coarse 
graining scale k much larger than the confinement scale guarantees that the complicated 
nonperturbative physics does not play a role yet. In this case, will look similar to the 
classical action, typically with a running gauge coupling evaluated at the scale k. (This 
does not hold for Green functions with much larger momenta p 2 ^> k 2 where the relevant 
IR cutoff is p, and the effective coupling is a s (p).) For lower k the resolution is smeared out 
and the detailed information of the short distance physics can be lost. (Again, this does not 
concern Green functions at high momenta.) On the other hand, the "observable volume" 
is increased and long distance aspects such as collective phenomena become visible. In a 
theory with a physical UV cutoff A we may associate T A with the classical action S since no 
fluctuations are effectively included. By definition, the effective average action equals the 
effective action for k = 0, Tq = T, since the infrared cutoff is absent. Thus r& interpolates 
between the classical action S and the effective action r as k is lowered from A to zero. 
The ability to follow the evolution k — > is equivalent to the ability to solve the quantum 
field theory. 

For a formal description we will consider in the first two sections a model with real scalar 
fields x a (the index a labeling internal degrees of freedom) in d Euclidean dimensions with 
classical action S[x\- We define the generating functional for the connected Green functions 
by 

W k [J] = In/ Vxexp^-Sklx] + J d d xJ a {x) X a {x)} (1.1) 
where we have added to the classical action an IR cutoff A^S* 

S k [x] = S[ X ]+AS k [ X ] (1.2) 
which is quadratic in the fields and conveniently formulated in momentum space 

ASk[x] = \S ^)- d R ^ 2 )Xa(-Q)x a (q) ■ (1-3) 
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Here J a are the usual scalar sources introduced to define generating functional and Rk(q 2 ) 
denotes an appropriately chosen (see below) IR cutoff function. We require that Rk(q 2 ) 
vanishes rapidly for q 2 ^> k 2 whereas for q 2 <C k 2 it behaves as Rk{q 2 ) — k 2 . This means 
that all Fourier components of \ a with momenta smaller than the IR cutoff k acquire an 
effective mass m e s ~ k and therefore decouple while the high momentum components of 
X a should not be affected by Rk- The classical fields 

#■ . <*•> = ^ta (i.4) 

now depend on k. In terms of Wk the effective average action is defined via a Legendre 
transformation 

r fc [$] = -W k [J] + J d d xJ a (x)$ a (x) - AS k m . (1.5) 

In order to define a reasonable coarse grained free energy we have subtracted in Eq.( |1.5| ) 
the infrared cutoff piece. This guarantees that the only difference between I\ and T is the 
effective IR cutoff in the fluctuations. Furthermore, this has the consequence that Tk does 
not need to be convex whereas a pure Legendre transform is always convex by definition. 
(The coarse grained free energy becomes convex || only for k — > 0.) This is very important 
for the description of phase transitions, in particular first order ones. One notes 

limi4(g 2 ) = => limr fc [$] = r[$] 

fc^o fc^o fl 6l 

lim R k (q 2 ) = oo lim T fc [$] = S[$\ 

k^A k^A 



for a convenient choice of Ru like 



R k (q 2 ) = Z^ k q 2 2/A 2 _ 2 /k 2 ■ (1.7) 



Here denotes the wave function renormalization to be specified below and we will 
often use for Rk the limit A — > oo 

Mf) = ■ (1-8) 

We note that the property Ta = 5* is not essential since the short distance laws may be 
parameterized by Ta as well as by S. In addition, for momentum scales much smaller than 
A universality implies that the precise form of Ta is irrelevant, up to the values of a few 
relevant renormalized couplings. 

A few properties of the effective average action are worth mentioning: 

1. All symmetries of the model which are respected by the IR cutoff AkS are auto- 
matically symmetries of Tk- In particular, this concerns translation and rotation 
invariance and one is not plagued by many of the problems encountered by a formu- 
lation of the block-spin action on a lattice. 
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2. In consequence, T k can be expanded in terms of invariants with respect to these 
symmetries with couplings depending on k. For the example of a scalar theory one 
may use a derivative expansion (p = $ a $ a /2) 

r k = J d d x \u k ( P ) + l -z^ k { P )d^ a d^ a + . . . j (1.9) 

and expand further in powers of p 

1_ 2 1_ 3 

£4(p) = i^k{p-po{k)) +Qlk(p-Po(k)) +••• 
^$,fc(p) = Z*,fc(po) + Z$ ifc (p ) (p - p ) + • • • (1-10) 

where p denotes the (fc-dependent) minimum of the effective average potential U k (p). 
We see that T k describes infinitely many running couplings. {Zq>^ in Eq.( |1.7|) can be 
identified with ^^(po)-) 

3. There is no problem incorporating chiral fermions since a chirally invariant cutoff R k 
can be formulated |K| [iT |. 

4. Gauge theories can be formulated along similar lines [[E^ II, 15, [15, [I7| even 
though A k S may not be gauge invariant. In this case the usual Ward identities 
receive corrections for which one can derive closed expressions [15] . These corrections 
vanish for k — > 0. 

5. The high momentum modes are very effectively integrated out because of the expo- 
nential decay of R k for q 2 ^> k 2 . Nevertheless, it is sometimes technically easier to 
use a cutoff without this fast decay property (e.g. R k ~ k 2 or R k ~ k 4 /q 2 ). In the 
latter case one has to be careful with possible remnants of an incomplete integration 
of the short distance modes. Also our cutoff does not introduce any non-analytical 
behavior as would be the case for a sharp cutoff 0. 

6. Despite a similar spirit and many analogies there remains also a conceptual difference 
to the Wilsonian effective action ■ The Wilsonian effective action describes a set 
of different actions (parameterized by A) for one and the same theory — the n-point 
functions are independent of A and have to be computed from by further func- 
tional integration. In contrast, T k describes the effective action for different theories 

- for any value of k the effective average action is related to the generating functional 
of a theory with a different action S k = S + A^S*. The n-point functions depend on 



k. The Wilsonian effective action does not generate the 1PI Green functions |]19 . 



7. Because of the incorporation of an infrared cutoff, T k is closely related to an effective 
action for averages of fields 0, where the average is taken over a volume ~ k d . 
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2 Exact renormalization group equation 



The dependence of T k on the coarse graining scale k is governed by an exact renormalization 
group equation (ERGE) [EO] 



d t r k [$} 



Tr 



dt.Ru 



[2.1) 



ln(fc/A) with some arbitrary momentum scale A, and the trace includes a mo- 

d d q 



Here t 

mentum integration as well as a summation over internal indices, Tr = / 



(2) 

second functional derivative T k denotes the exact inverse propagator 



(2tt)< 



,(2) 



ab 



Ea- The 



(2.2) 



The flow equation Eq. (|2.1|) can be derived from Eq.( |1.5| ) in a straightforward way using 



dtT k \t = -d t W k \j-d t A k S[$) 



and 



5 2 W k 



-Tr{d t R k ((®$) 
lTr{d t R k W®} 



5 2 (T k + A k S) 



<d>> ($))} 



(2.3) 



S 2 W k 



5Ja(-q)5J b (q') 5$ b (-q')5$ c (q") 



5J a (-q)5J b (q') 



ftacSqq" 



It has the form of a renormalization group improved one-loop expression 
one-loop formula for T k reads 

r fc [$] = £[d>] + ~Trln (S^m + R k 



(2.4) 
Indeed, the 

(2.5) 



with the second functional derivative of the classical action, similar to Eq. (|2.2j ). (Re- 
member that is the field dependent classical inverse propagator. Its first and second 
derivative with respect to the fields describe the classical three- and four-point vertices, 
respectively.) Taking a t-derivative of Eq.( |2.5| ) gives a one-loop flow equation very similar 
to Eq.Q with rf } replaced by S^. It may seem surprising, but it is nevertheless true, 



that the renormalization group improvement — > T k promotes the one-loop flow equa- 
tion to an exact nonperturbative flow equation which includes the effects from all loops 
as well as all contributions which are non-analytical in the couplings like instantons, etc.! 
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For practical computations it is actually often quite convenient to write the flow equation 
Eq. ( |2.1|) as a formal derivative of a renormalization group improved one-loop expression 

d t r k =~Trd t ln(r£ ] +R k ) (2.6) 

with dt acting only on R k and not on T k , i.e. d t = (dR k /dt) (d/dR k ). Flow equations 
for n-point functions follow from appropriate functional derivatives of Eq. ( [2.1] ) or Eq. ( |2.6| ) 
with respect to the fields. For their derivation it is sufficient to evaluate the corresponding 
one-loop expressions (with the vertices and propagators derived from T k ) and then to take 
a formal ^-derivative. (If the one-loop expression is finite or properly regularized the dt~ 
derivative can be taken after the evaluation of the trace.) This permits the use of (one-loop) 
Feynman diagrams and standard perturbative techniques in many circumstances. Most 
importantly, it establishes a very direct connection between the solution of flow-equations 
and perturbation theory. If one uses on the right hand side of Eq.( |2.1| ) a truncation for 
which the propagator and vertices appearing in Y k are replaced by the ones derived from 
the classical action, but with running fc-dependent couplings, and then expands the result 
to lowest non-trivial order in the coupling constants one recovers standard renormalization 
group improved one-loop perturbation theory. The formal solution of the flow equation 
can also be employed for the development of a systematically resummed perturbation 
theory pi . 



For a choice of the cutoff function similar to Eq.( |l.&| ) the momentum integral contained 
in the trace on the right hand side of the flow equation is both infrared and ultraviolet 
finite. Infrared finiteness arises through the presence of the infrared regulator ~ R k . We 

(2,) 

note that all eigenvalues of the matrix T k + R k must be positive semi-definite. The proof 
follows from the observation that the functional T k + A k S is convex since it is obtained from 
W k by a Legendre transform. On the other hand, ultraviolet finiteness is related to the fast 
decay of dtRk for q 2 3> k 2 . This expresses the fact that only a narrow range of fluctuations 
with q 2 ~ k 2 contributes effectively if the infrared cutoff k is lowered by a small amount. If 
for some other choice of R k the right hand side of the flow equation would not remain UV 
finite this would indicate that the high momentum modes have not yet been integrated out 
completely in the computation of r*.. Since the flow equation is manifestly finite this can 
be used to define a regularization scheme. The "ERGE-scheme" is specified by the flow 
equation, the choice of R k and the "initial condition" T\. This is particularly important 
for gauge theories where other regularizations in four dimensions and in the presence of 
chiral fermions are difficult to construct. For gauge theories Ta has to obey appropriately 
modified Ward identities. In the context of perturbation theory a first proposal how to 
regularize gauge theories by use of flow equations can be found in fl3| , |14}] . We note that 
in contrast to previous versions of exact renormalization group equations there is no need 
in the present formulation to construct an ultraviolet momentum cutoff — a task known 
to be very difficult in non-Abelian gauge theories. 

Despite the conceptual differences between the Wilsonian effective action and the 
effective average action T k the exact flow equations describing the A-dependence of S 1 ^ and 
the fc-dependence of T k are simply related. Polchinski's continuum version of the Wilsonian 
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flow equation 0] can be transformed into Eq.( |2.1| ) by means of a Legendre transform and 
a suitable variable redefinition |22| . 



Even though intuitively simple, the replacement of the (RG-improved) classical propa- 
gator by the full propagator turns the solution of the flow equation Eq. (|2.1| ) into a difficult 
mathematical problem: The evolution equation is a functional differential equation. Once 
T k is expanded in terms of invariants (e.g. Eqs. (|1.9|) , (|1.10|) ) this is equivalent to a coupled 
system of non-linear partial differential equations for infinitely many couplings. General 
methods for the solution of functional differential equations are not developed very far. 
They are mainly restricted to iterative procedures that can be applied once some small 
expansion parameter is identified. This covers usual perturbation theory in the case of 
a small coupling, the 1/iV-expansion or expansions in the dimensionality 4 — d or 2 — d. 
It may also be extended to less familiar expansions like a derivative expansion which is 
related in critical three dimensional scalar theories to a small anomalous dimension. In 
the absence of a clearly identified small parameter one nevertheless needs to truncate the 
most general form of 1^ in order to reduce the infinite system of coupled differential equa- 
tions to a (numerically) manageable size. This truncation is crucial. It is at this level 
that approximations have to be made and, as for all nonperturbative analytical methods, 
they are often not easy to control. The challenge for nonperturbative systems like low 
momentum QCD is to find flow equations which (a) incorporate all the relevant dynamics 
such that neglected effects make only small changes, and (b) remain of manageable size. 
The difficulty with the first task is a reliable estimate of the error. For the second task the 
main limitation is a practical restriction for numerical solutions of differential equations 
to functions depending only on a small number of variables. The existence of an exact 
functional differential flow equation is a very useful starting point and guide for this task. 
At this point the precise form of the exact flow equation is quite important. Furthermore, 
it can be used for systematic expansions through enlargement of the truncation and for 
an error estimate in this way. Nevertheless, this is not all. Usually, physical insight into a 
model is necessary to device a useful nonperturbative truncation! 

So far, two complementary approaches to nonperturbative truncations have been ex- 
plored: an expansion of the effective Lagrangian in powers of derivatives (p = |$ a $ a ) 

r fc [$] = J d d x \u k ( P ) + l -z^)d^ a d^ a + -Y^Md^p + o(<9 4 ) j (2.7) 

or one in powers of the fields 

00 -1 „ / n \ 

r*[*] = £ - / II Pfo) - $ o] r k n \ Xl , ...,x n ). (2.8) 

n=0 n - J \j=0 ) 



If one chooses $ as the fc-dependent VEV of $, the series Eq.( [2.8| ) starts effectively at 
n — 2. The flow equations for the 1PI n-point functions 1^ are obtained by functional 
differentiation of Eq.( J2.1[ ). Such flow equations have been discussed earlier from a somewhat 
different viewpoint 0. They can also be interpreted as a differential form of Schwinger- 
Dyson equations . 
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The formation of mesonic bound states, which typically appear as poles in the (Min- 
kowskian) four-quark Green function, is most efficiently described by expansions like 
Eq.( [2.8| ). This is also the form needed to compute the nonperturbative momentum de- 
pendence of the gluon propagator and the heavy quark potential [[U|, [17|, On the other 
hand, a parameterization of r& as in Eq.( j2.7Q seems particularly suited for the study of 
phase transitions. The evolution equation for the average potential Uk follows by evaluat- 
ing Eq. (|2.7| ) for constant $. In the limit where the ^-dependence of Z$^ is neglected and 
y$,fc = one finds ]7) for the 0(iV)-symmetric scalar model 

B TJ ( ) = I [ ^_^± ( jV ~ 1 + \ ^ (O Q) 

* k[P) 2 J (2nY dt \Z^ h q* + R k + U' k + Z^q 2 + R„{q) + U' k + 2pU^J 1 " ] 

with U' k = ^p 1 , etc. One observes the appearance of p-dependent mass terms in the effective 
propagators of the right hand side of Eq . (|2~9| ) . Once ^$ = — d t lnZ^ t k is determined 
in terms of the couplings parameterizing Uk this is a partial differential equation for a 
function Uk depending on two variables k and p which can be solved numerically |24|, |25], 
|2"7| . (The Wilson-Fisher fixed point relevant for a second order phase transition id = 3) 
corresponds to a scaling solution |28| |2^] where d t U k = 0.) A suitable truncation of a flow 



equation of the type Eq.( [2.7|) will play a central role in the description of chiral symmetry 
breaking below. 

It should be mentioned at this point that the weakest point in the ERGE approach seems 
to be a reliable estimate of the truncation error in a nonperturbative context. This problem 
is common to all known analytical approaches to nonperturbative phenomena and appears 
often even within systematic (perturbative) expansions. One may hope that the existence of 
an exact flow equation could also be of some help for error estimates. An obvious possibility 
to test a given truncation is its enlargement to include more variables — for example, 
going one step higher in a derivative expansion. This is similar to computing higher orders 
in perturbation theory and limited by practical considerations. As an alternative, one 
may employ different truncations of comparable size — for instance, by using different 
definitions of retained couplings. A comparison of the results can give a reasonable picture 
of the uncertainty if the used set of truncations is wide enough. In this context we should 
also note the dependence of the results on the choice of the cutoff function Rk(q)- Of 
course, for k — >• the physics should not depend on a particular choice of Rk and, in fact, 
it does not for full solutions of Eq.([2.1|). Different choices of Rk just correspond to different 
trajectories in the space of effective average actions along which the unique IR limit r[$] 
is reached. Once approximations are used to solve the ERGE Eq. (|2.1| ), however, not only 
the trajectory but also its end point will depend on the precise definition of the function 
Rk- This is very similar to the renormalization scheme dependence usually encountered in 
perturbative computations of Green functions. One may use this scheme dependence as a 
tool to study the robustness of a given approximation scheme. 

Before applying a new nonperturbative method to a complicated theory like QCD it 
should be tested for simpler models. A good criterion for the capability of the ERGE to 
deal with nonperturbative phenomena concerns the critical behavior in three dimensional 
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scalar theories. In a first step the well known results of other methods for the critical 
exponents have been reproduced within a few percent accuracy [^|. The ability of the 
method to produce new results has been demonstrated by the computation of the critical 
equation of state for Ising and Heisenberg models which has been verified by lattice 
simulations |3(| . This has been extended to first order transitions in matrix models [|27| or 
for the Abelian Higgs model relevant for superconductors |3T| , p6| . Analytical investigations 
of the high temperature phase transitions in d — 4 scalar theories (O(iV)-models) have 
correctly described the second order nature of the transition in contrast to earlier 
attempts within high temperature perturbation theory. 

For an extension of the flow equations to Abelian and non-Abelian gauge theories we 
refer the reader to |L2|, [13], [14], [ER |lj| |17|, |i~8| . The other necessary piece for a description 



of low-energy QCD, namely the transition from fundamental (quark and gluon) degrees 
of freedom to composite (meson) fields within the framework of the ERGE can be found 
in [33|. We will describe the most important aspects of this formalism for mesons below. 



3 Chiral symmetry breaking in QCD 

The strong interaction dynamics of quarks and gluons at short distances or high energies 
is successfully described by quantum chromodynamics (QCD). One of its most striking 
features is asymptotic freedom j34J which makes perturbative calculations reliable in the 
high energy regime. On the other hand, at scales around a few hundred MeV confinement 
sets in. As a consequence, the low-energy degrees of freedom in strong interaction physics 
are mesons, baryons and glueballs rather than quarks and gluons. When constructing 
effective models for these IR degrees of freedom one usually relies on the symmetries of 
QCD as a guiding principle, since a direct derivation of such models from QCD is still 
missing. The most important symmetry of QCD, its local color SU(3) invariance, is of not 
much help here, since the IR spectrum appears to be color neutral. When dealing with 
bound states involving heavy quarks the so called "heavy quark symmetry" may be invoked 
to obtain approximate symmetry relations between IR observables ||35|| . We will rather 
focus here on the light scalar and pseudoscalar meson spectrum and therefore consider 
QCD with only the light quark flavors u, d and s. To a good approximation the masses 
of these three flavors can be considered as small in comparison with other typical strong 
interaction scales. One may therefore consider the chiral limit of QCD (vanishing current 
quark masses) in which the classical QCD Lagrangian does not couple left- and right- 
handed quarks. It therefore exhibits a global chiral invariance under Ul(N) x Ur(N) = 
SUl(N) x SUr(N) x Uy (1) x £/a(1) where N denotes the number of massless quarks (N = 2 
or 3) which transform as 

Ub$r\ U R eU R (N) 

U L ^ L - U L eU L {N). (3.1) 



yR = — 2 — ^ 

/ - 1 + ^ , 
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Even for vanishing quark masses only the diagonal SUv{N) vector-like subgroup can be 
observed in the hadron spectrum ("eightfold way"). The symmetry SUl(N) x SUr(N) 
must therefore be spontaneously broken to SUv{N) 

SU L (N) x SU R (N) — > SU L+R (N) = SU V (N) . (3.2) 

Chiral symmetry breaking is one of the most prominent features of strong interaction 



dynamics and phenomenologically well established [pqj , though a rigorous derivation of 
this phenomenon starting from first principles is still missing. In particular, the chiral 
symmetry breaking Eq.(|3.2|) predicts for N = 3 the existence of eight light parity-odd 



(pseudo-)Goldstone bosons: tt°, tt^, K°, K° , K ± and 77. Their comparably small masses 
are a consequence of the explicit chiral symmetry breaking due to small but non-vanishing 
current quark masses. The axial Abelian subgroup U^(l) = Ul_r(1) is broken in the quan- 
tum theory by an anomaly of the axial-vector current. This breaking proceeds without the 
occurrence of a Goldstone boson [[T/[]. Finally, the Uy(l) = Ul+r(1) subgroup corresponds 



to baryon number conservation. 

The light pseudoscalar and scalar mesons are thought of as color neutral quark-anti- 
quark bound states $ afe ~ ip L ip R , a,b = 1, . . . , N, which therefore transform under chiral 
rotations Eq. ( |3.1| ) as 

$ — ► u R m\ . (3.3) 



Hence, the chiral symmetry breaking pattern Eq.( [3.2|) is realized if the meson potential 
develops a VEV 

'X b ) = a S ab ; Zfo^O. (3.4) 



One of the most crucial and yet unsolved problems of strong interaction dynamics is to 
derive an effective field theory for the mesonic degrees of freedom directly from QCD which 
exhibits this behavior. 



4 A semi— quantitative picture 

Before turning to a quantitative description of chiral symmetry breaking using flow equa- 
tions it is perhaps helpful to give a brief overview of the relevant scales which appear in 
relation to this phenomenon and the physical degrees of freedom associated to them. Some 
of this will be explained in more detail in the remainder of these lectures whereas other 
parts are rather well established features of strong interaction physics. 

At scales above approximately 1.5 GeV, the relevant degrees of freedom of strong in- 
teractions are quarks and gluons and their dynamics appears to be well described by per- 
turbative QCD. At somewhat lower energies this changes dramatically. Quark and gluon 
bound states form and confinement sets in. Concentrating on the physics of scalar and 
pseudoscalar mesons [] there are three important momentum scales which appear to be 
rather well separated: 

1 One may assume that all other bound states are integrated out. We will comment on this issue below. 
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• The compositeness scale k$ at which mesonic ipip bound states form because of the 
increasing strength of the strong interaction. It will turn out to be somewhere in the 



• The chiral symmetry breaking scale k x $B at which the chiral condensate (ip tjj a ) or 



$ ) assumes a non-vanishing value, therefore breaking chiral symmetry according 



to Eq.(|372|). This scale is found to be around (400 — 500) MeV. For k below k xSB the 
quarks acquire constituent masses M q ~ 350 MeV due to their Yukawa coupling to 



• The confinement scale Aqcd — 200 MeV which corresponds to the Landau pole in 
the perturbative evolution of the strong coupling constant a s . In our context, this is 
the scale where possible deviations of the effective quark propagator from its classical 
form and multi-quark interactions not included in the meson physics may become 
very important. 

For scales k in the range k x sB ^ k ^ the most relevant degrees of freedom are mesons 
and quarks. Typically, the dynamics in this range is dominated by the strong Yukawa 
coupling h between quarks and mesons: h 2 /(A7i) ^> a s . One may therefore assume that 
the dominant QCD effects are included in the meson physics and consider a simple model 
of quarks and mesons only. As one evolves to scales below k x s B the Yukawa coupling 
decreases whereas a s increases. Of course, getting closer to Aqcd it is no longer justified 
to neglect the QCD effects which go beyond the dynamics of effective meson degrees of 
freedom. On the other hand, the final IR value of the Yukawa coupling h is fixed by the 
typical values of constituent quark masses M q ~ 350 MeV to be h 2 /(47r) ~ 4.5. One may 
therefore speculate that the domination of the Yukawa interaction persists down to scales 
k — M q at which the quarks decouple from the evolution of the mesonic degrees of freedom 
altogether due to their mass. Of course, details of the gluonic interactions are expected to 
be crucial for an understanding of quark and gluon confinement. Strong interaction effects 
may dramatically change the momentum dependence of the quark n-point functions for k 
around Aqcd- Yet, as long as one is only interested in the dynamics of the mesons one is led 
to expect that these effects are quantitatively no too important. Because of the effective 
decoupling of the quarks and therefore the whole colored sector the details of confinement 
have only little influence on the mesonic flow equations for k ^ Aq C d- We conclude that 
there are good prospects that the meson physics can be described by an effective action 
for mesons and quarks for k < The main part of the work presented here is concerned 
with this effective quark meson model. 

In order to obtain this effective action at the compositeness scale kq> from short distance 
QCD two steps have to be carried out. In a first step one computes at the scale k p ~ 
1.5 GeV an effective action involving only quarks. This step integrates out the gluon 
degrees of freedom in a "quenched approximation" . More precisely, one solves a truncated 
flow equation for QCD with quark and gluon degrees of freedom in presence of an effective 
infrared cutoff k p ~ 1.5 GeV in the quark propagators. The exact flow equation to be 



range (600 - 700) MeV. 





the chiral condensate Eq. (|3.4j ). 
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Figure 1: Heavy quark potential for different values of a*. The solid line corresponds to 
a* = 1.5 whereas the dashed lines are for a* = 1 (lower dashed line) and a* = 2 (upper 
dashed line). Also shown are the two loop potential and the Richardson fit as dotted lines. 



used for this purpose is obtained by lowering the infrared cutoff Rk for the gluons to zero 
while keeping the one for the quarks fixed. Subsequently, the gluons are eliminated by 
solving the field equations for the gluon fields as functionals of the quarks. This will result 
in a non-trivial momentum dependence of the quark propagator and effective non-local 
four and higher quark interactions. Because of the infrared cutoff k p the resulting effective 
action for the quarks resembles closely the one for heavy quarks (at least for Euclidean 
momenta). The dominant effect is the appearance of an effective quark potential (similar 
to the one for the charm quark) which describes the effective four-quark interactions. For 
the effective quark action at k p we only retain this four-quark interaction in addition to 
the two-point function, while neglecting n-point functions involving six and more quarks. 

For typical momenta larger than Aqcd a reliable computation of the effective quark 
action should be possible by using in the quark-gluon flow equation a relatively simple 
truncation. The result |1| for the Fourier transform of the potential V(q 2 ) is shown 
in figure ffl. (See ref for the relation between the four-quark interaction and the 
heavy quark potential which involves a rescaling in dependence on k p .) A measure for the 
truncation uncertainties is the parameter a* which corresponds to the value to which an 
appropriately defined running strong coupling a s evolves for k — > 0. We see that these 
uncertainties affect principally the low-g 2 region. For high values of q 2 the potential is 
very close to the perturbative two-loop potential whereas for intermediate q 2 relevant for 
quarkonium spectra it is quite close to phenomenologically acceptable potentials (e.g. the 
Richardson potential p8|). The inverse quark propagator is found in this computation to 
remain very well approximated by the simple classical momentum dependence ^ . 
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In the second step one has to lower the infrared cutoff in the effective non-local quark 
model in order to extrapolate from k p to This task can be carried out by means of the 
exact flow equation for quarks only, starting at k p with an initial value Tk p [ip] as obtained 
after integrating out the gluons. For fermions the trace in Eq.( |2.1|) has to be replaced by 
a supertrace in order to account for the minus sign related to Grassmann variables pi . A 



first investigation in this direction |33[ has used a truncation with a chirally invariant four 
quark interaction whose most general momentum dependence was retained 



r, 



+ 




V;(p)^, fc (p) i>5 ab + m a \p) lb + im a \p) MP) 



X ^k\pi,P2,P3,pd{\fi' a (-pM(P2)} [^b(P^j(-P3 



S(-Pi)75^i(p 2 )] {4(P*h^(-Pz)} } • (4-1) 



Here i,j run from one to N c which is the number of quark colors. The indices a, b label the 
different light quark flavors and run from 1 to N. The matrices m and m are hermitian and 
m + 277175 forms therefore the most general quark mass matrix. (Our chiral conventions |10 



where the hermitean part of the mass matrix is multiplied by 75 may be somewhat unusual 
but they are quite convenient for Euclidean calculations.) The ansatz Eq. (|4.1|) does not 
correspond to the most general chirally invariant four-quark interaction. It neglects similar 
interactions in the p-meson and pomeron channels which are also obtained from a Fierz 
transformation of the heavy quark potential With V(q 2 ) the heavy quark potential in 
a Fourier representation, the initial value at k p = 1.5 GeV was taken as (Z^ k — Z^ k(p 2 = 

-K)) 



y f t ) (Pl,P2,P3,P4)Z 



2 

ip,k p 



-V((pi-p 3 ) 



+ 



8ttA 



(Pi-Ps) 2 ((pi-Ps) 2 Y 



(4.2) 



This corresponds to an approximation by a one gluon exchange term ~ <y s (k p ) and a 
string tension A ~ 0.18 GeV 2 and is in reasonable agreement with the form computed 
recently |1| from the solution of flow equations (see figjl]). In the simplified ansatz Eq. ( |4.2j ) 
the string tension introduces a second scale in addition to k p and it becomes clear that 
the incorporation of gluon fluctuations is a crucial ingredient for the emergence of mesonic 



bound states. For a more precise treatment |L8[ of the four-quark interaction at the scale 
k$ this second scale is set by the running of ot s or Aqcd- 

The evolution equation for the function A^. for k < k p can be derived from the fermionic 

version of Eq. (pTT|) and the truncation Eq. (|4.1|). Since A^. depends on six independent 
momentum invariants it is a partial differential equation for a function depending on seven 
variables and has to be solved numerically 



The "initial value" Eq. 



corresponds to 



the t-channel exchange of a "dressed" colored gluonic state and it is by far not clear that 
the evolution of A^ will lead at lower scales to a momentum dependence representing the 
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exchange of colorless mesonic bound states. Yet, at the compositeness scale 



fc$~ 630MeV (4.3) 

one finds an approximate factorization 

a£ } (pi,P2,P3,P4) = g(pi,P2)G(s)g(p3,pi) + ... (4.4) 

which indicates the formation of mesonic bound states. Here g{p\,P2) denotes the ampu- 
tated Bethe-Salpeter wave function and G(s) is the mesonic bound state propagator dis- 
playing a pole-like structure in the s-channel if it is continued to negative s = (pi +P2) 2 - 
The dots indicate the part of Ajj. which does not factorize and which will be neglected in 
the following. In the limit where the momentum dependence of g and G is neglected we re- 



cover the four-quark interaction of the Nambu-Jona-Lasinio model |3j|, ^(J . It is therefore 
not surprising that our description of the dynamics for k < k<p will parallel certain aspects 
of other investigations of this model, even though we are not bound to the approximations 
used typically in such studies (large-iV c expansion, perturbative renormalization group, 
etc.). 

It is clear that for scales k ^ kq, a description of strong interaction physics in terms of 
quark fields alone would be rather inefficient. Finding physically reasonable truncations 
of the effective average action should be much easier once composite fields for the mesons 
are introduced. The exact renormalization group equation can indeed be supplemented by 
an exact formalism for the introduction of composite field variables or, more generally, a 
change of variables [33]. For our purpose, this amounts in practice to inserting at the scale 
kq, the identities 

1 = const J VaA 

x exp j- tr (a A - K A G -m A - ] 6) i (a A - GK A -m A - GO) } 

1 = const J Va H (4.5) 

x exp j - tr (a\j - K ] H G - m ] H - C (5)t G) i (a H - GK H - m H - GC (5) ) j 

into the functional integral which formally defines the quark effective average action. Here 
we have used the shorthand notation A^GB = J -^ ![ A* a {q)G ah {q)B} ) {q)^ and Ka,h are 
sources for the collective fields <ja,h which correspond in turn to the anti-hermitian and 
hermitian parts of the meson field $. They are associated to the fermion bilinear operators 
0[ip], 0^[ijj] whose Fourier components read 

j^y9(-p,p + q)*P a (p)Mp + q) 

{5)a b (q) = - J -^jl9(-P,P + q)i> a {p)l^b{p + q) ■ (4.6) 
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The choice of g{—p,p + q) as the bound state wave function renormalization and of G(q) as 
its propagator guarantees that the four-quark interaction contained in Eq.( |4.5|) cancels the 
dominant factorizing part of the QCD-induced non-local four-quark interaction Eqs. (|4.1|) , 
( |4.4|) . In addition, one may choose 

m T Hab = m ab (0)g- 1 (0,0)Z i;M X0) 

™ T Aab = m ab (0)g-\0,0)Z^ ki X0) (4.7) 

such that the explicit quark mass term cancels out for q = 0. The remaining quark bilinear 
is ~ m(q) —m(0)Z^^(0)g(—q, q)/[Z^p }k ^(q)g(0, 0)]. It vanishes for zero momentum and will 
be neglected in the following. Without loss of generality we can take m real and diagonal 
and m = 0. 

In consequence, we have replaced at the scale the effective quark action Eq. (|4.1| ) 
with Eq.( |4.4j ) by an effective quark meson action given by 

?k = r fe -i|rf 4 xtr($tj + y$) 

T k = Jd 4 xU k ($,&) (4.8) 
d 4 q 



+ 



(27T) 



Z^ k (q)q 2 tr [$ f (g)$(g)] + Z^q^^q^q^q) 



At the scale /c$ the inverse scalar propagator is related to G(q) in Eq.( |4.4j ) by 

G- 1 (q 2 ) = 2ml+2Z< s>M (q)q 2 . (4.9) 

This fixes the term in U k<s> which is quadratic in <3> to be positive, U k<s> = fn\ tr*^^ + . . .. 
The higher order terms in U ki> cannot be determined in the approximation Eq.( [4.1| ) since 
they correspond to terms involving six or more quark fields. The initial value of the Yukawa 
coupling corresponds to the "quark wave function in the meson" in Eq.([4.4|), i.e. 



hk*(-q,q-p) = g(-q,q-p) (4.10) 

which can be normalized with h kiS> (0,0) = g(0,0) = 1. We observe that the explicit chiral 
symmetry breaking from non-vanishing current quark masses appears now in the form of 
a meson source term with 

J = 2m^Z v , iA , $ (0)5' _1 (0,0) (m ab + im ab ) = 2Z^ M ml <s> diag(m u ,m d ,m s ) . (4.11) 

This induces a non-vanishing ($) and an effective quark mass M q through the Yukawa 
coupling. We note that the current quark mass m q and the constituent quark mass M q ~ 
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hk ($) are identical at the scale kq>. (By solving the field equation for $ as a functional of 
ip, ip (with Uk = m|tr$^$) one recovers from Eq. ( [4.8|) the effective quark action Eq . (|4 . 1|) . 
For a generalization beyond the approximation of a four-quark interaction or a quadratic 
potential see ref [01 •) Spontaneous chiral symmetry breaking can be described in this 
language by a non-vanishing ($) in the limit j — > 0. Because of spontaneous chiral 
symmetry breaking the constituent quark mass M q can differ from zero even for m q = 0. It 
is a nice feature of our formalism that it provides for a unified description of the concepts 
of the current and the constituent quark masses. As long as the effective average potential 
Uk has a unique minimum at $ = there is simply no difference between the two. The 
running of the current quark mass in the pure quark model should be equivalent in the 
quark meson language to the running of hk ($) (k). (A verification of this property would 
actually provide a good check for the truncation errors.) Nevertheless, the formalism is 
now adapted to account for the quark mass contribution from chiral symmetry breaking 
since the absolute minimum of Uk may be far from the perturbative one for small k. We 
also note that Tk in Eq.( |4.8|) is chirally invariant. The explicit chiral symmetry breaking 
in Tk appears only in the form of a linear source term which is independent of k and does 
not affect the flow of I\. The flow equation for Tk therefore respects chiral symmetry even 
in the presence of quark masses. This leads to a considerable simplification. 

At the scale fc$ the propagator G and the wave function g(—q, q—p) should be optimized 
for a most complete elimination of terms quartic in the quark fields. In the present context 
we will, however, neglect the momentum dependence of Z^^, Z<$,^ and hk- The mass 
was found in |33| for the simple truncation Eq.( |4.1|) with Z^^ = 1, m = m = to be 



— 120 MeV. In view of the possible large truncation errors we will take this only as 
an order of magnitude estimate. Below we will consider the range = (45 — 120) MeV 
for which chiral symmetry breaking can be obtained in a two flavor model. Furthermore, 
we will assume, as usually done in large-iV c computations within the NJL-model, that 
Z<b,k<s, = Z^^ifl = 0) < 1. The quark wave function renormalization = Z^^(q = 0) is 
set to one at the scale kq> for convenience. For k < kq> we will therefore study an effective 
action for quarks and mesons in the truncation 

r fc = J d^z^j^r + z^kti [d^d^] + u k {<s>, $ f ) 

+ ^(^^-^(^V 6 } (4-12) 
with compositeness conditions 

U k ^{<5>,&) = ml $ tr$t$- -z7 fc# (det $ + det $ f ) 

1- / + \2 N - 1- / + 1 + \ 2 
+ (tr$ f $) +^^A 2 , fc$ tr ^$ t $-— tr^^J +... 

= -J— ~ (45 MeV) 2 - (120 MeV) 2 
2G(0) 
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hk 9 — Z^, )ki> — 1 

Z* M < 1. (4.13) 
As a consequence, the initial value of the renormalized Yukawa coupling which is given 

— -. 1/2 

by hk 9 = hk^Z^^Z^J^ is large! Note that we have included in the potential an explicit 
Ua{X) breaking term ~ Vk which mimics the effect of the chiral anomaly of QCD to 
leading order in an expansion of the effective potential in powers of Because of the 
infrared stability of the evolution of Tk which will be discussed below the precise form of 
the potential, i.e. the values of the quartic couplings and so on, will turn out to be 
unimportant. 

We have refrained here for simplicity from considering four quark operators with vector 
and pseudo-vector spin structure. Their inclusion is straightforward and would lead to 
vector and pseudo-vector mesons in the effective action Eq.( 4.12 ). We will concentrate 



first on two flavors and consider only the two limiting cases Vj. = and — > oo. We 
also omit first the explicit quark masses and study the chiral limit j — 0. Because of the 
positive mass term m| $ one has at the scale kq, a vanishing expectation value (<3>) = 
(for j = 0). There is no spontaneous chiral symmetry breaking at the compositeness scale. 
This means that the mesonic bound states at kq, and somewhat below are not directly 
connected to chiral symmetry breaking. 

The question remains how chiral symmetry is broken. We will try to answer it by 
following the evolution of the effective potential JJ\. from k$, to lower scales using the 
exact renormalization group method outlined above with the compositeness conditions 
Eq.( [4.13D defining the initial values. In this context it is important that the formalism 



for composite fields |33[ also induces an infrared cutoff in the meson propagator. The 
flow equations are therefore exactly of the form Eq.( [2.1| ) (except for the supertrace), with 
quarks and mesons treated on an equal footing. In fact, one would expect that the large 
renormalized Yukawa coupling will rapidly drive the scalar mass term to negative values 
as the IR cutoff k is lowered ||10|| . This will then finally lead to a potential minimum away 
from the origin at some scale k x sB < kq> such that ($) ^ 0. The ultimate goal of such a 
procedure, besides from establishing the onset of chiral symmetry breaking, would be to 
extract phenomenological quantities, like or meson masses, which can be computed in 
a straightforward manner from in the IR limit k — > 0. 

At first sight, a reliable computation of r fc ^ seems a very difficult task. Without a 
truncation r& is described by an infinite number of parameters (couplings, wave function 
renormalizations, etc.) as can be seen if 1^ is expanded in powers of fields and derivatives. 
For instance, the pseudoscalar and scalar meson masses are obtained as the poles of the 
exact propagator, lim fc ^ r^ 2 ' ) (g)|$ = ^), which receives formally contributions from terms 
in Tfc with arbitrarily high powers of derivatives and the expectation value a . Realistic 
nonperturbative truncations of Tk which reduce the problem to a manageable size are 
crucial. We will argue in the following that there may be a twofold solution to this problem: 

• Due to an IR fixed point structure of the flow equations in the symmetric regime, 



2 The anomaly term in the fermionic effective average action has been computed in |4f 
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i.e. for k x sB < k < the values of many parameters of r& for k — > will be 
approximately independent of their initial values at the compositeness scale For 
small enough Zq>^ only a few relevant parameters (ml , V^) need to be computed 
accurately from QCD. They can alternatively be determined from phenomenology. 

One can show that physical observables like meson masses, decay constants, etc., can 
be expanded in powers of the quark masses within the linear meson model [^2], [43], [44 . 



This is similar to the way it is usually done in chiral perturbation theory [B5|. To a 
given finite order of this expansion only a finite number of terms of a simultaneous 
expansion of in powers of derivatives and $ are required if the expansion point is 
chosen properly. 

In combination, these two results open the possibility for a perhaps unexpected degree of 
predictive power within the linear meson model. 

We wish to stress, though, that a perturbative treatment of the model at hand, e.g., 
using perturbative RG techniques, cannot be expected to yield reliable results. The renor- 
malized Yukawa coupling is expected to be large at the scale Even the IR value of h 
is still relatively big 

lim/i=^~7.5 (4.14) 

and h increases with increasing k. The dynamics of the linear meson model is therefore 
clearly nonperturbative for all scales k < k$. 



5 Flow equations for the linear quark meson model 

We will next turn to the ERGE analysis of the linear meson model which was introduced in 
the last section?]. In a first approach we will attack the problem at hand by truncating 
in such a way that it contains all perturbatively relevant and marginal operators, i.e. those 
with canonical dimensions. This has the advantage that the flow of a small number of cou- 
plings permits quantitative insight into the relevant mechanisms, e.g., of chiral symmetry 
breaking. More quantitative precision will be obtained once we generalize our truncation 
(see below) for the 0(4)-model by allowing for the most general effective average potential 
Uk- The effective potential is a function of only four SUl(N) x SUr(N) invariants for 
N = 3: 

p = tr$ f $ 

N /^t* 1 \ 2 
r 2 = _tr(*t* ,j 

t 3 = tr(W__p) 

f = det$ + det$ f . (5.1) 

3 For a study of chiral symmetry breaking in QED using related exact renormalization group techniques 
see ref [Q. 
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The invariant T3 is only independent for N > 3. For iV = 2 it can be eliminated by a suitable 
combination of r 2 and p. The additional Ua(1) breaking invariant u = z(det $ — det $') 
is CP violating and may therefore appear only quadratically in Uk- It is straightforward 
to see that uj 2 is expressible in terms of the invariants Eq.( |5.1|) . We may expand Uk as a 
function of these invariants around its minimum, i.e. p 
t~2 = T3 = where 



po = Nal £ = £0 = 2of and 



0-Q 



(5.2) 



The expansion coefficients are the fc-dependent couplings of the model. In our first version 
we only keep couplings of canonical dimension d c < 4. This yields in the chirally symmetric 
regime, i.e., for /c x sb < k < k® where Wq = 



U k 



m lp 



*V 2 W-1 T 



2,fc r 2 



2^ 



whereas in the SSB regime for < /c x sb we have 



£4 — 2^ 1 ' k 



P 



2 JV-1- 1_ 
H - — A 2 ,a,.t 2 + -z^fc 



■iV-2 
0,k 



(5.3) 



(5.4) 



Before continuing to compute the nonperturbative beta functions for these couplings it 
is worthwhile to pause here and emphasize that naively (perturbatively) irrelevant opera- 
tors can by no means always be neglected. The most prominent example for this is QCD 
itself. It is the very assumption of our treatment of chiral symmetry breaking (substanti- 



ated by the results of p3| ) that the momentum dependence of the coupling constants of 
some six-dimensional quark operators {ipip) 2 develop poles in the s-channel indicating the 
formation of mesonic bound states. On the other hand, it is quite natural to assume that 
$ 6 or $ 8 operators are not really necessary to understand the properties of the potential in 
a neighborhood around its minimum. Yet, truncating higher dimensional operators does 
not imply the assumption that the corresponding coupling constants are small. In fact, this 
could only be expected as long as the relevant and marginal couplings are small as well. 
What is required, though, is that their influence on the evolution of those couplings kept 
in the truncation, for instance, the set of equations Eq.( |5.6| ) below, is small. A comparison 
with the results of the later sections for the full potential Uk (i.e., including arbitrarily 
many couplings in the formal expansion) will provide a good check for the validity of 
this assumption. In this context it is perhaps also interesting to note that the truncation 
Eq. (4.12) includes the known one-loop beta functions of a small coupling expansion as well 



as the leading order result of the large- iV c expansion of the Ul(N) x Ur(N) model [46 



This should provide at least some minimal control over this truncation, even though we 
believe that our results are significantly more accurate. 

Inserting the truncation Eqs. fl4.12D , ( [5. 3D , (|5.4|) into Eq. (|2.1j ) reduces this functional 
differential equation for infinitely many variables to a finite set of ordinary differential 
equations. This yields, in particular, the beta functions for the couplings Ai^, \2 ) k, 
and m 2 k or 0o,fc- Details of the calculation can be found in WJ§. We will refrain here from 
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presenting the full set of flow equations but rather illustrate the main results with a few 
examples. Defining dimensionless renormalized VEV and coupling constants 



K 
,2 



Z^ jk Nal k k 2 

Ir 

Aj = \i : kZ$ 2 k ; i = 1, 2 
— — 

v = V k Z^l k " 4 

one finds, e.g., for the spontaneous symmetry breaking (SSB) regime and z7 fc = 



d 
dt 



Jt 1 



dt 



-(2 + ??$)/? + 



1 



N 2 li(0;riz) + 3£(2Ai«; 77*) 



+ (iV 2 -l) 
= 2r]$Ai + 



1 + 



167T 2 

A 2 



Ai 



h 2 



16tt 2 



AT AjZf (0; r/$) + 9A?Z£(2Ai«; 77$) 



+ (iV 2 - 1) [A x + A 2 ] 2 4(A 2 k; ^) - Aj^hH^^^K; ifc) j 

1 f 7V 2 9 
= 2r/ $ A 2 + — _A 2 4(0;^) + ^(iV 2 - 4)A 2 Z 4 (A 2 *;; 77*) 



,4;(F)4, 



(5.5) 



(5.6) 



±h 2 
dt 



-N 2 X 2 2 lf x (0, A 2 k; I/*) + 3[A 2 + 4Ax]A 2 / 4 1 (2A 1 /t, A 2 k; 77$) 



[d - 4 + 27/0 + r/„] /i 2 - -^ d /i 4 { AT 2 zg B)d (l^ 2 , 0; ^, 77*) 



- (iy 2 -l)C )d (-^ 2 , K A 2 ;^,r/ $ ) 

Here 77$ = — In 77^, = — J- In Z^ jk are the meson and quark anomalous dimensions, re- 
spectively f47fl . The symbols 1%, I* n2 and l^ d denote bosonic and fermionic mass threshold 
functions, respectively, which are defined in ||47|| . They describe the decoupling of massive 
modes and provide an important nonperturbative ingredient. For instance, the bosonic 
threshold functions 



pi, .... .. , " ' !//2/< d 



d d q 1 dR k 



(2n) d Z^ k dt [P(g2) + A;2 



w 



in+1 



(5.7) 



involve the inverse average propagator P(q 2 ) = q 2 + Z^ l k R k (q 2 ) where the infrared cutoff is 
manifest. These functions decrease ~ w — for w ^> 1. Since typically w = M 2 /k 2 with 
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M a mass of the model, the main effect of the threshold functions is to cut off fluctuations 
of particles with masses M 2 ^> k 2 . Once the scale k is changed below a certain mass 
threshold, the corresponding particle no longer contributes to the evolution and decouples 
smoothly. 

Within our truncation the beta functions Eq.( |5.6| ) for the dimensionless couplings look 
almost the same as in one-loop perturbation theory. There are, however, two major new 
ingredients which are crucial for our approach: First, there is a new equation for the 
running of the mass term in the symmetric regime or for the running of the potential 
minimum in the regime with spontaneous symmetry breaking. This equation is related 
to the quadratic divergence of the mass term in perturbation theory and does not appear 
in the Callan-Symanzik [[f8j or Coleman-Weinberg j49j treatment of the renormalization 
group. Obviously, these equations are the key for a study of the onset of spontaneous 
chiral symmetry breaking as k is lowered from k$ to zero. Second, the most important 
nonperturbative ingredient in the flow equations for the dimensionless Yukawa and scalar 
couplings is the appearance of effective mass threshold functions like Eq. (|5.7|) which account 
for the decoupling of modes with masses larger than k. Their form is different for the 
symmetric regime (massless fermions, massive scalars) or the regime with spontaneous 
symmetry breaking (massive fermions, massless Goldstone bosons). Without the inclusion 
of the threshold effects the running of the couplings would never stop and no sensible 
limit k — > could be obtained because of unphysical infrared divergences. The threshold 
functions are not arbitrary but have to be computed carefully. The mass terms appearing 
in these functions involve the dimensionless couplings. Expanding the threshold functions 
in powers of the mass terms (or the dimensionless couplings) makes their non-perturbative 
content immediately visible. It is these threshold functions which will make it possible 
below to use one-loop type formulae for the necessarily nonperturbative computation of 
the critical behavior of the (effectively) three-dimensional 0(4)-symmetric scalar model. 



6 The chiral anomaly and the O (4)— model 

We have seen how the mass threshold functions in the flow equations describe the de- 
coupling of heavy modes from the evolution of as the IR cutoff k is lowered. In the 
chiral limit with two massless quark flavors (N = 2) the pions are the massless Gold- 
stone bosons. Below, once we will consider the flow of the full effective average potential 
without resorting to a quartic truncation, we will also consider the more general case of 
non-vanishing current quark masses. For the time being, the effect of the physical pion 
mass of m n ~ 140 MeV, or equivalently of the two small but non-vanishing current quark 
masses, can easily be mimicked by stopping the flow of at k = m w by hand. This 
situation changes significantly once the strange quark is included. Now the r\ and the four 
K mesons appear as additional massless Goldstone modes in the spectrum. They would 
artificially drive the running of at scales ^ k ^ 500 MeV where they should already 
be decoupled because of their physical masses. It is therefore advisable to focus on the 
two flavor case N = 2 as long as the chiral limit of vanishing current quark masses is 
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considered. 

It is straightforward to obtain an estimate of the (renormalized) coupling v which 
parameterizes the explicit £/a(1) breaking due to the chiral anomaly. From Eq. (|5.4j) we 
find 

m 2, = jua^- 2 ~ 1 GeV (6.1) 

which translates for N = 2 for k — > into 

u(k -»• 0) ~ 1 GeV . (6.2) 

This suggests that v — ► oo can be considered as a realistic limit. An important simplifi- 
cation occurs for N = 2 and f — ► oo, related to the fact that for N = 2 the chiral group 
SUl{2) x SUr{2) is (locally) isomorphic to 0(4). Thus, the complex (2,2) representation 
$ of SUl{2) x SUr{2) may be decomposed into two vector representations, (a, 7i k ) and 
(r]',a k ) of 0(4): 

$ = i(<T-i7/) + I( a * + Mr fe )T fc . (6.3) 

For v — > oo the masses of the 7/ and the a fc diverge and these particles decouple. We are 
then left with the original 0(4) symmetric linear a-model of Gell-Mann and Levy |]5D| 
coupled to quarks. The flow equations of this model have been derived previously |H], [7| 
for the truncation of the effective action used here. For mere comparison we also consider 
the opposite limit v — *■ 0. Here the rf meson becomes an additional Goldstone boson in 
the chiral limit which suffers from the same problem as the K and the r\ in the case N = 3. 
Hence, we may compare the results for two different approximate limits of the effects of 
the chiral anomaly: 

• the 0(4) model corresponding to N = 2 and v — > oo 

• the Ul(2) x Ur{2) model corresponding to N = 2 and v = 0. 

For the reasons given above we expect the first situation to be closer to reality. In this 
case we may imagine that the fluctuations of the kaons, r], 7/ and the scalar mesons (as 
well as vector and pseudovector mesons) have been integrated out in order to obtain the 
initial values of — in close analogy to the integration of the gluons for the effective 
quark action [-0] discussed above. We will keep the initial values of the couplings , 
Ai fe.j, and Z$ >kit as free parameters. Our results should be quantitatively accurate to the 
extent to which the local polynomial truncation is a good approximation. 

7 Infrared stability 

Eq. (|5.6|) and the corresponding set of flow equations for the symmetric regime constitute 
a coupled system of ordinary differential equations which can be integrated numerically. 
Similar equations can be computed for the 0(4) model where N = 2 and the coupling 
A2 is absent. We have neglected in a first step the dependence of all threshold functions 
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Figure 2: Evolution of the renormalized mass m in the symmetric regime (dashed line) and 
the vacuum expectation value cxo = Z^Wq^ of the scalar field in the SSB regime (solid 
line) as functions of k for the f/z,(2) x Ur(2) model. Initial values at k = kq> are = 0, 
hj = 300 and m fc$ = 63 MeV. 



appearing in the flow equations on the anomalous dimensions. This dependence will be 
taken into account below once we abandon the quartic potential approximation. The most 
important result is that chiral symmetry breaking indeed occurs for a wide range of initial 
values of the parameters including the presumably realistic case of large renormalized 
Yukawa coupling and a bare mass of order 100 MeV. A typical evolution of the 
renormalized meson mass m with k is plotted in fig. p[ Driven by the strong Yukawa 
coupling, m decreases rapidly and goes through zero at a scale /c x sb not far below k$. Here 
the system enters the SSB regime and a non-vanishing (renormalized) VEV o"o for the 
meson field $ develops. The evolution of o"o with k turns out to be reasonably stable already 
before scales k ~ where the evolution is stopped. We take this result as an indication 
that our truncation of the effective action r\ leads at least qualitatively to a satisfactory 
description of chiral symmetry breaking. The reason for the relative stability of the IR 
behavior of the VEV (and all other couplings) is that the quarks acquire a constituent 
mass M q = hao ~ 350 MeV in the SSB regime. As a consequence they decouple once 
k becomes smaller than M q and the evolution is then dominantly driven by the massless 
Goldstone bosons. This is also important in view of potential confinement effects expected 
to become important for the quark n-point functions for k around Aqcd — 200 MeV. Since 
confinement is not explicitly included in our truncation of T^, one might be worried that 
such effects could spoil our results completely. Yet, as discussed in some more detail above, 
only the colored quarks should feel confinement and they are no longer important for the 
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evolution of the meson couplings for k around 200 MeV. One might therefore hope that 
a precise treatment of confinement is not crucial for this approach to chiral symmetry 
breaking. 

Most importantly, one finds that the system of flow equations exhibits a partial IR fixed 
point in the symmetric phase. As already pointed out one expects Z$ t k$ to be rather small. 
In turn, one may assume that, at least for the initial range of running in the symmetric 
regime the mass parameter m 2 ~ is large. This means, in particular, that all threshold 
functions with arguments ~ m 2 may be neglected in this regime. As a consequence, the 
flow equations simplify considerably. We find, for instance, for the Ul{2) x Ur{2) model 
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(7.1) 



This system possesses an attractive IR fixed point for the quartic scalar self interactions 

2 

N ' 



Ai* 



X A 

-A 2 * 



(7.2) 



Furthermore it is exactly soluble f47fl . Because of the strong Yukawa coupling the quartic 
couplings Ai and A 2 generally approach their fixed point values rapidly, long before the 
systems enters the broken phase (m — > 0) and the approximation of large m breaks down. 
In addition, for large initial values h 2 the Yukawa coupling at the scale k x sB where m 
vanishes (or becomes small) only depends on the initial value m^. Hence, the system 
is approximately independent in the IR upon the initial values of Ai, A 2 and h 2 , the 
only "relevant" parameter being m^. (Once quark masses and a proper treatment of 
the chiral anomaly are included for iV = 3 one expects that m q and v are additional 
relevant parameters. Their values may be fixed by using the masses of tt, K and rf 
as phenomenological input.) In other words, the effective action looses almost all its 
"memory" in the far IR of where in the UV it came from. This feature of the flow equations 
leads to a perhaps surprising degree of predictive power. In addition, also the dependence of 
f n = 2o"o on mfc # is not very strong for a large range in , as shown in fig. ||. The relevant 
parameter m k<i can be fixed by using the constituent quark mass M q = hcr ~ 350 MeV as 
a phenomenological input. One obtains for the 0(4)-model 
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The resulting value for the decay constant is 

f v = (91 - 100) MeV 



(7.3) 
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Figure 3: The pion decay constant f w as a function of m^ $ /A;| for = 630 MeV and initial 
values \j = and hf = 300 (solid line) as well as = 10 4 (dashed line). 



for h 2 j = 10 4 — 300. It is striking that this comes close to the real value f w = 92.4 MeV but 
we expect that the uncertainty in the determination of the compositeness scale and the 
truncation errors exceed the influence of the variation of hj. We have furthermore used 
this result for an estimate of the chiral condensate: 



™2 



-m 



f n Z-y* A~-(195MeVf 



(7.5) 



where the factor A ~ 1.7 accounts for the change of the normalization scale of ypipj from 
kq, to the commonly used value 1 GeV. Our value is in reasonable agreement with results 
from sum rules |5T|. This is non-trivial since not only and f w enter but also the IR 
value Z$ t k=o- Integrating Eq. (|7.1| ) for 77$ one finds 



— ln T 



(7.6) 



Thus Z$ t k w iU indeed be practically independent of its initial value Z^^ already after 
some running as long as Z<$,^ is small compared to 0.01. 

The alert reader may have noticed that the beta functions Eq.( [7.1| ) correspond exactly 
to those obtained in the one-quark-loop approximation or, in other words, to the leading 
order in the large-iV c expansion for the Nambu-Jona-Lasinio model [I5|. The fixed point 
Eq.( [7.2| ) is then nothing but the large-iV c boundary condition on the evolution of Ai and 
A2 in this model. Yet, we wish to stress that nowhere we have made the assumption 
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that N c is a large number. On the contrary, the physical value N c = 3 suggests that the 
large-iV c expansion should a priori only be trusted on a quantitatively rather crude level. 
The reason why we expect Eq. ([7.1|) to nevertheless give rather reliable results is based on 
the fact that for small Zq>^ all (renormalized) meson masses are much larger than the 
scale k for the initial part of the running. This implies that the mesons are effectively 
decoupled and their contribution to the beta functions is negligible leading to the one- 
quark-loop approximation. Yet, already after some relatively short period of running the 
renormalized meson masses approach zero and our approximation of neglecting mesonic 
threshold functions breaks down. Hence, the one-quark-loop approximation is reasonable 
only for scales close to k$ but is bound to become inaccurate around k x sB and in the 
SSB regime. What is important in our context is not the numerical value of the partial 
fixed points Eq. ( |7.2| ) but rather their mere existence and the presence of a large coupling 
h 2 driving the Aj fast towards them. This is enough for the IR values of all couplings to 
become almost independent of the initial values Xij. Similar features of IR stability are 
expected if the truncation is enlarged, for instance, to a more general form of the effective 
potential as will be discussed in the next section. 



8 Flow equation for the scalar potential 

In this and the remaining sections we consider the 0(4)-symmetric quark meson model 
without truncating its effective average potential to a polynomial form |52| . A comparison 
with the results of the last sections will give us a feeling for the size of the errors induced 
by the quartic truncation used so far. Furthermore, such an approach is well suited for 
a study of the chiral phase transition close to which the form of the potential deviates 
substantially from a polynomial. 

It is convenient to work with dimensionless and renormalized variables therefore elimi- 
nating all explicit fc-dependence. With (t = ln(fc/fc$)) 

u(t,p) = k- d U k {p) , p = Z^ k k 2 ~ d p (8.1) 

and using Eq. (|4.12|) as a first truncation of the effective average action T k one obtains the 
flow equation in arbitrary dimensions d 

d 

—u = —du + (d — 2 + 77$) pu' 

at (g_ 2 ) 
+ 2v d |3/ V; V*) + itW + 2pu"; m ) - 2l +1 N c l { F)d ^ph 2 ; % )} . 

Here v^ 1 = 2 d+1 n d ^ 2 T(d/2) and primes denote derivatives with respect to p. We will 
always use in the following for the number of quark colors iV c = 3. Eq. (|8.2|) is a partial 
differential equation for the effective potential u(t, p) which has to be supplemented by the 
flow equation for the Yukawa coupling and expressions for the anomalous dimensions 77$, 
77^,. The definition of the threshold functions l d , can be found in prTfl . 
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The dimensionless renormalized expectation value k = 2k 2 d Z<§,j.o\ k , with o"o,fc the 
fc-dependent VEV of $, may be computed for each k directly from the condition 

u'(t, «) = -J=k~^Zlf = e g (8.3) 
V2k 



where 52 



J = ( 8 - 4 ) 

and rh = (m u + m^)/2 denotes the average light current quark mass normalized at 
Note that k = in the symmetric regime for vanishing source term. Eq. ( |8.3| ) allows us to 
follow the flow of k according to 



with A = u"(t, k). We define the Yukawa coupling for p = k and its flow equation reads [47] 
-h 2 = (d - 4 + 2?fy + 77$) /i 2 - 2v d h 4 ^3l[ F 1 B)d (^h 2 K,e g ; 71^,7]$) 



d_ 

dt' 



A F l B)d (^h 2 K, 6 g + 2\K] ffy, //$) 



8.6) 



Similarly, the the scalar and quark anomalous dimensions are infered from 



V<s> = -^\nZ$ )k = A^\AK\ 2 m d 22 (e g ,6 g + 2\K]ri<b) 



+ 2iN c h 2 m[ F) \\h 2 ^) 



5.7) 



+ m[ F 2 B)d {^h 2 K,e g + 2\K-,r] i ,,r]^) 



which is a linear set of equations for the anomalous dimensions. The definitions of the 
threshold functions l^^f, nT>n iri2 , m& d and mj^^ are again specified in 



The flow equations (8.2), (|8.5| ) — (|8.7|) , constitute a coupled system of ordinary and 
partial differential equations which can be integrated numerically |)3|, 54jl . Here we take 
the effective current quark mass dependence of h, Z$ t k and Z^ k into account by stopping 
the evolution according to Eqs. (|8.6| ), (|8.7]) , evaluated for the chiral limit, below the pion 
mass tUtt. 

Similarly to the case of the quartic truncation of the effective average potential described 
in the preceding section one finds for d = 4 that chiral symmetry breaking indeed occurs 
for a wide range of initial values of the parameters. These include the presumably realistic 
case of large renormalized Yukawa coupling and a bare mass of order 100 MeV. 
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Most importantly, the approximate partial IR fixed point behavior encountered for the 
quartic potential approximation before, carries over to the truncation of which maintains 



the full effective average potential [52]. To see this explicitly we study the flow equations 



?D subject to the condition <C 1. For the relevant range of p both 

u'(t,p) and u'(t,p) + 2pu"(t,p) are then much larger than ph 2 (t) and we may therefore 
neglect in the flow equations all scalar contributions with threshold functions involving 
these large masses. This yields the simplified equations (d = 4, v^ 1 = 32ii 2 ) 
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Of 

dt 



— Au + (2 + 77$) pv! 

8tt 2 ' 
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jOkiWfloh 2 ) 



(8.8) 
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Again, this approximation is only valid for the initial range of running below before 
the (dimensionless) renormalized scalar mass squared u'(t,p = 0) approaches zero near the 
chiral symmetry breaking scale. The system Eq.(|8.8|) is exactly soluble. We find 



h 2 (t) = Z^{t) 



hj 



1 



u{t,p) 



e- At uj(e 2t p h2it) 



N c ft 



hj 



^ /o dre-H^\\h\t)pe 
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™2r\ 



(The integration over r on the right hand side of the solution for u can be carried out by 
first exchanging it with the one over momentum implicit in the definition of the threshold 
function Iq F ^.) Here Ui(p) = w(0,/5) denotes the effective average potential at the compos- 
iteness scale and hj is the initial value of h 2 at /c$, i.e. for t = 0. For simplicity we will use 
an expansion of the initial value effective potential uj(p) in powers of p around p = 
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(8.10) 



even though this is not essential for the forthcoming reasoning. Expanding also f )4 m 
Eq. (|8.9|) in powers of its argument one finds for n > 2 
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For decreasing t — > — oo the initial values uf^ become rapidly unimportant and vS n ^ /h 2n 
approaches a fixed point. For n = 2, i.e., for the quartic coupling, one finds 



u( 2 >(t,0) 
h 2 (t) 



1 



4 2) (o) 



M<Lh 2 f 



(8.12) 
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leading to the fixed point value (u^/h 2 )* = 1 already encountered in Eq. ([7.21). As a 
consequence of this fixed point behavior the system looses all its "memory" on the initial 
values Uj 2 ^ at the compositeness scale This typically happens before the approxima- 
tion u'(t, p),u'(t, p) + 2pu"(t,p) ^> ph 2 {t) breaks down and the solution Eq. (|S.9|) becomes 
invalid. Furthermore, the attraction to partial infrared fixed points continues also for the 
range of k where the scalar fluctuations cannot be neglected anymore. As for the quartic 
truncation the initial value for the bare dimensionless mass parameter 



u'j(0) fn 



•2 



h 2 k\ 



(8.13) 



is never negligible. (In fact, using the values for mjL and k<$, computed previously |33| for 
a pure quark effective action as described above one obtains rri\ it /k% ~ 0.036.) For large 
hj (and dropping the constant piece w/(0)) the solution Eq. QS.9D therefore behaves with 
growing \t\ as 

- ~ Wt ~ (8-14) 

hj 27l 2 Jo 2 

In other words, for hi — > oo the IR behavior of the linear quark meson model will depend 
(in addition to the value of the compositeness scale kq> and the quark mass rh) only on one 
parameter, m| . We have numerically verified this feature by starting with different values 

for uf\o). Indeed, the differences in the physical observables were found to be small. This 
IR stability of the flow equations leads to a perhaps surprising degree of predictive power: 
Not only the scalar wave function renormalization but even the full effective potential 
U(p) = lim fc ^ Uk(p) is (approximately) fixed for Z$ k<s> 1 once is known! For 
definiteness we will perform our numerical analysis of the full system of flow equations 
(|S.2|), (|8.5|) — (|8.7| ) with the idealized initial value ui(p) = W/(0)p in the limit hj — > oo. It 
should be stressed, though, that deviations from this idealization will lead only to small 
numerical deviations in the IR behavior of the linear quark meson model as long as the 



condition <S 1 holds, say for [42| hj ^ 15. 

With this knowledge at hand we may now fix the remaining three parameters of our 
model, k$, m? and rh by using f T = 92.4 MeV, the pion mass M n = 135 MeV and the 
constituent quark mass M q as phenomenological input. This approach differs from that 
of the preceding sections where we took an earlier determination of as input for the 
computation of f v . It will be better suited for precision estimates of the high temperature 
behavior in the following sections. Because of the uncertainty regarding the precise value of 
M q we give in tab. [I] the results for several values of M q . The first line of tab. p] corresponds 
to the choice of M q and A/ = u"(k) which we will use for the forthcoming analysis of the 
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-(180) 3 


80.5 


339 
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0.0225 


64.4 


19.5 


13.7 


-(174) 3 
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Table 1: The table shows the dependence on the constituent quark mass M q of the input parameters 
fn^/k% and j as well as some of our "predictions". The phenomenological input used here besides 
M q is /tt = 92.4 MeV, m T = 135 MeV. The first line corresponds to the values for M q and A/ used in the 
remainder of this work. The other three lines demonstrate the insensitivity of our results with respect to 
the precise values of these parameters. 



model at finite temperature. As argued analytically above the dependence on the value of 
Xj is weak for large enough hi as demonstrated numerically by the second line. Moreover, 
we notice that our results, and in particular the value of j, are rather insensitive with 
respect to the precise value of M q . It is remarkable that the values for fc$ and are not 



very different from those computed ]33| from four-quark interactions as described above. 
As compared to the analysis of the preceding sections the present truncation of is of 
a higher level of accuracy: We now consider an arbitrary form of the effective average 
potential instead of a polynomial approximation and we have included the pieces in the 
threshold functions which are proportional to the anomalous dimensions. It is encouraging 
that the results are rather robust with respect to these improvements. 

Once the parameters k$, mj^ and rh are fixed there are a number of "predictions" of 
the linear meson model which can be compared with the results obtained by other methods 
or direct experimental observation. First of all one may compute the value of rh at a scale 
of 1 GeV which is suitable for comparison with results obtained from chiral perturbation 
theory pH] and sum rules [^1[]. For this purpose one has to account for the running of this 
quantity with the normalization scale from as given in tab. |]to the commonly used value 
of 1 GeV: m(l GeV) = A rh(k$). A reasonable estimate of the factor A is obtained from 
the three loop running of rh in the MS scheme For M q ~ 300 MeV corresponding to 
the first two lines in tab. [I] its value is A ~ 1.3. The results for rh(l GeV) are in acceptable 
agreement with recent results from other methods l)R |51j even though they tend to be 



somewhat larger. Closely related to this is the value of the chiral condensate 

,-1/2 



(1 GeV) = -Ami f« Z *,k=o - 2m 



(8.15) 



These results are quite non-trivial since not only f n and m\ enter but also the computed 
IR value Z^>^ =Q . We emphasize in this context that there may be substantial corrections 
both in the extrapolation from /c$ to 1 GeV and because of the neglected influence of the 
strange quark which may be important near /c$. These uncertainties have only little effect 
on the physics at lower scales as relevant for our analysis of the temperature effects. Only 
the value of j which is fixed by m n enters here. 
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A further more qualitative test concerns the mass of the sigma resonance or radial mode 
in the limit k — > whose renormalized mass squared is given by 



ml = Z^ 2 -^ + 4Aa 2 . (8.16) 

From our numerical analysis we obtain Xk=o — 20 which translates into m a ~ 430 MeV. 
One should note, though, that this result is presumably not very accurate as we have 
employed in this work the approximation of using the Goldstone boson wave function 
renormalization constant also for the radial mode. Furthermore, the explicit chiral sym- 
metry breaking contribution to m 2 a is certainly underestimated as long as the strange quark 
is neglected. In any case, we observe that the sigma meson is significantly heavier than the 
pions. This is a crucial consistency check for the linear quark meson model. A low sigma 
mass would be in conflict with the numerous successes of chiral perturbation theory p6 | 



which requires the decoupling of all modes other than the Goldstone bosons in the IR-limit 
of QCD. The decoupling of the sigma meson is, of course, equivalent to the limit A — > oo 
which formally describes the transition from the linear to the non-linear sigma model and 
which appears to be reasonably well realized by the large IR-values of A obtained in our 
analysis. We also note that the issue of the sigma mass is closely connected to the value 
of f(°\ the value of f n in the chiral limit m = also given in tab. [1]. To lowest order in 
{fir — )//tt or > equivalently, in m one has 



Z][ ml m 



A larger value of m a would therefore reduce the difference between ffi and f w . 

In fig. ^ we show the dependence of the pion mass and decay constant on the average 
current quark mass m. These curves depend very little on the values of the initial param- 
eters as demonstrated in tab. |TJ by f^. We observe a relatively large difference of 12 MeV 
between the pion decay constants at rh — m p h ys and m — 0. According to Eq. (|8.17|) this 



difference is related to the mass of the sigma particle and will be modified in the three 
flavor case. We will later find that the critical temperature T c for the second order phase 
transition in the chiral limit is almost independent of the initial conditions. The values 
of and T c essentially determine the non-universal amplitudes in the critical scaling 
region (see below). In summary, we find that the behavior of our model for small k is quite 
robust as far as uncertainties in the initial conditions at the scale are concerned. We will 
see that the difference of observables between non-vanishing and vanishing temperature is 
entirely determined by the flow of couplings in the range < k ^ 3T. 



9 Chiral phase transition of two flavor QCD 

Strong interactions in thermal equilibrium at high temperature T — as realized in early 
stages of the evolution of the Universe — differ in important aspects from the well tested 
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vacuum or zero temperature properties. A phase transition at some critical temperature 
T c or a relatively sharp crossover may separate the high and low temperature physics |)6] . 
Many experimental activities at heavy ion colliders |57| search for signs of such a transition. 
It was realized early that the transition should be closely related to a qualitative change 
in the chiral condensate according to the general observation that spontaneous symmetry 
breaking tends to be absent in a high temperature situation. A series of stimulating 
contributions []5B|, [5I| |5t| pointed out that for sufficiently small up and down quark masses, 
m u and m^, and for a sufficiently large mass of the strange quark, m s , the chiral transition 
is expected to belong to the universality class of the 0(4) Heisenberg model. This means 
that near the critical temperature only the pions and the sigma particle play a role for the 
behavior of condensates and long distance correlation functions. It was suggested |59|, |60] 
that a large correlation length may be responsible for important fluctuations or lead to 

El 



This was even related |59], [6(J to the spectacular 



a disoriented chiral condensate 
"Centauro events" [B2| observed in cosmic rays. The question how small m u and would 
have to be in order to see a large correlation length near T c and if this scenario could be 
realized for realistic values of the current quark masses remained, however, unanswered. 
The reason was the missing link between the universal behavior near T c and zero current 
quark mass on one hand and the known physical properties at T = for realistic quark 
masses on the other hand. 

It is the purpose of the remaining sections of these lectures to provide this link []52 



We present here the equation of state for two flavor QCD within an effective quark me- 
son model. The equation of state expresses the chiral condensate O^'ip) as a function 
of temperature and the average current quark mass rh = (m u + m^)/2. This connects 
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Figure 5: The plot shows the chiral condensate (ipip) as a function of temperature T. Lines (a), (6), (c), 
(d) correspond at zero temperature to m n — 0, 45 MeV, 135 MeV, 230 MeV, respectively. For each pair of 
curves the lower one represents the full T-dependence of (ipip) whereas the upper one shows for comparison 
the universal scaling form of the equation of state for the 0(4) Heisenberg model. The critical temperature 
for zero quark mass is T c = 100.7 MeV. The chiral condensate is normalized at a scale fc$ ~ 620 MeV. 

explicitly the universal critical behavior for T — > T c and rh — > with the temperature 
dependence for a realistic value m phys . Since our discussion covers the whole tempera- 
ture range < T£,1.7T C we can fix m p h ys such that the (zero temperature) pion mass is 
= 135 MeV. The condensate (ipip^ plays here the role of an order parameter, fig. |5| 

shows our results for (rf>ip\ (T, rh): Curve (a) gives the temperature dependence of {^4^ 
in the chiral limit rh = 0. Here the lower curve is the full result for arbitrary T whereas 
the upper curve corresponds to the universal scaling form of the equation of state for the 
0(4) Heisenberg model. We see perfect agreement of both curves for T sufficiently close to 
T c = 100.7 MeV. This demonstrates the capability of our method to cover the critical be- 
havior and, in particular, to reproduce the critical exponents of the 0(4)-model. We have 
determined (see below) the universal critical equation of state as well as the non-universal 
amplitudes. This provides the full functional dependence of (ipip^ (T, rh) for small T — T c 
and rh. The curves (b), (c) and (d) are for non-vanishing values of the average current 
quark mass rh. Curve (c) corresponds to m p h ys or, equivalently, m^(T = 0) = 135 MeV. 
One observes a crossover in the range T = (1.2 — 1.5)T C . The 0(4) universal equation 
of state (upper curve) gives a reasonable approximation in this temperature range. The 
transition turns out to be much less dramatic than for rh = 0. We have also plotted in 
curve (b) the results for comparably small quark masses ~ 1 MeV, i.e. rh = m p h ys /10, for 
which the T = value of equals 45 MeV. The crossover is considerably sharper but a 
substantial deviation from the chiral limit remains even for such small values of rh. In order 
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m,7r 
MeV 





45 


135 


230 


T m 

J pc 

MeV 


100.7 


~ 110 


~ 130 


~ 150 


T W 

pc 

MeV 


100.7 


113 


128 





Table 2: The table shows the critical and "pseudocritical" temperatures for various values of the zero 
temperature pion mass. Here Tp\) is defined as the inflection point of (iptp) (T) whereas Tpc is the location 
of the maximum of the sigma correlation length. 



to facilitate comparison with lattice simulations which are typically performed for larger 
values of m n we also present results for m n (T = 0) = 230 MeV in curve (d). One may 
define a "pseudocritical temperature" T pc associated to the smooth crossover as the inflec- 
tion point of (ipip^ (T) as often done in lattice simulations. Our results for this definition 

of T pc are denoted by T^) and are presented in tab. |2| for the four different values of rh or, 
equivalently, m n (T = 0). The value for the pseudocritical temperature for = 230 MeV 
compares well with the lattice results for two flavor QCD (cf. the discussion below). One 
should mention, though, that a determination of T pc according to this definition is subject 
to sizeable numerical uncertainties for large pion masses as the curve in fig. [5] is almost 
linear around the inflection point for quite a large temperature range. A difficult point in 
lattice simulations with large quark masses is the extrapolation to realistic values of m n or 
even to the chiral limit. Our results may serve here as an analytic guide. The overall pic- 
ture shows the approximate validity of the 0(4) scaling behavior over a large temperature 
interval in the vicinity of and above T c once the (non-universal) amplitudes are properly 
computed. 

A second important result of our investigations is the temperature dependence of the 
space-like pion correlation length m~ 1 (T). (We will often call m n (T) the temperature 
dependent pion mass since it coincides with the physical pion mass for T = 0.) The plot 
for m n (T) in fig. [] again shows the second order phase transition in the chiral limit rh — 0. 
For T < T c the pions are massless Goldstone bosons whereas for T > T c they form together 
with the sigma particle a degenerate vector of 0(4) whose mass increases as a function 
of temperature. For rh = the behavior for small positive T — T c is characterized by 
the critical exponent i.e. (T) = (i + Y l T c ({T -T c )/T c y and we obtain v = 0.787, 
£ + = 0.270. For rh > we find that m 7r (T) remains almost constant for T ^ T c with only 
a very slight dip for T near T c /2. For T > T c the correlation length decreases rapidly and 
for T ^> T c the precise value of rh becomes irrelevant. We see that the universal critical 
behavior near T c is quite smoothly connected to T = 0. The full functional dependence 
of m n (T,m) allows us to compute the overall size of the pion correlation length near the 
critical temperature and we find m^Tpc) ~ 1.7m 7r (0) for the realistic value m p h ys . This 
correlation length is even smaller than the vacuum (T = 0) one and shows no indication for 
strong fluctuations of pions with long wavelength. It would be interesting to see if a decrease 
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Figure 6: The plot shows m w as a function of temperature T for three different values of the average light 
current quark mass to. The solid line corresponds to the realistic value to = TO p hy S whereas the dotted line 
represents the situation without explicit chiral symmetry breaking, i.e., rh = 0. The intermediate, dashed 
line assumes to = TO p h ys /10. 



of the pion correlation length at and above T c is experimentally observable. It should be 
emphasized, however, that a tricritical behavior with a massless excitation remains possible 
for three flavors. This would not be characterized by the universal behavior of the 0(4)- 
model. We also point out that the present investigation for the two flavor case does not 
take into account a speculative "effective restoration" of the axial £/a(1) symmetry at high 
temperature |58L B3j. We will comment on these issues in the last section. In the next 



sections we will describe the formalism which leads to these results. 



10 Thermal equilibrium and dimensional reduction 

The extension of flow equations to thermal equilibrium situations at non-vanishing temper- 
ature T is straightforward |32| . In the Euclidean formalism non-zero temperature results 
in (anti-)periodic boundary conditions for (fermionic) bosonic fields in the Euclidean time 
direction with periodicity |64| 1/T. This leads to the replacement 

/(g J /(, 2 )-^E/(0/(* 2 (O + f 2 ) do,) 

in the trace in Eq.( [2.1|) when represented as a momentum integration, with a discrete 
spectrum for the zero component 

m _ / 2lnT for bosons , , 

qo[ > \ (2/ + 1)ttT for fermions. 1 } 
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Hence, for T > a four-dimensional QFT can be interpreted as a three-dimensional model 
with each bosonic or fermionic degree of freedom now coming in an infinite number of copies 
labeled by I G 2Z (Matsubara modes). Each mode acquires an additional temperature 
dependent effective mass term q^Q). In a high temperature situation where all massive 
Matsubara modes decouple from the dynamics of the system one therefore expects to 
observe an effective three-dimensional theory with the bosonic zero modes as the only 
relevant degrees of freedom. In other words, if the characteristic length scale associated 
with the physical system is much larger than the inverse temperature the compactified 
Euclidean "time" dimension cannot be resolved anymore. This phenomenon is known as 



'dimensional reduction" 65 



The formalism of the effective average action automatically provides the tools for a 
smooth decoupling of the massive Matsubara modes as the scale k is lowered from k ^> T 
to k <C T. It therefore allows us to directly link the low-T, four-dimensional QFT to the 
effective three-dimensional high-T theory. The replacement Eq.( |10.1| ) in Eq.( |2.1| ) mani- 
fests itself in the flow equations ( |8.2| ), ( p.5|) — ( |B.7|) only through a change to T-dependent 
threshold functions. For instance, the dimensionless functions l d ,(w;rj^) defined in Eq.( |5.7| ) 
are replaced by 



fc" ' 4 d (2 7 r) d - 1 \Z* tk dt J [P(q^) + k 2 w} n+1 

(10-3) 

where q 2 = q 2 + q 2 and go — 2ttIT. A list of the various temperature dependent threshold 
functions appearing in the flow equations can be found fl5"!2fl . There we also discuss some 
subtleties regarding the definition of the Yukawa coupling and the anomalous dimensions 
for T ^ 0. In the limit k ^> T the sum over Matsubara modes approaches the integration 
over a continuous range of go and we recover the zero temperature threshold function 
ln{w; 77$). In the opposite limit k <T the massive Matsubara modes (/ 7^ 0) are suppressed 
and we expect to find a d— 1 dimensional behavior of l d . In fact, one obtains from Eq.( |10.3| ) 

l d n (w,T/k; V<s> ) ~ li(w;rfr) for T < k , 

T flO 4) 

l d n (w,T/k-r)z) ~ -— lir^w;^) for T > k . 
k v d 



For our choice of the infrared cutoff function R^, Eq. (|1.8|) , the temperature dependent 
Matsubara modes in l d (w,T/k;r)$) are exponentially suppressed for T k whereas the 
behavior is more complicated for other threshold functions appearing in the flow equations 
( |S.2| ), ( |8.5f ) — (|8.7| ). Nevertheless, all bosonic threshold functions are proportional to T/k 
for T ^> k whereas those with fermionic contributions vanish in this limit[J This behav- 
ior is demonstrated in fig. ^ where we have plotted the quotients l\ (w, T/k)/l\(w) and 
lf ] \w,T/k)/l[ F) \w) of bosonic and fermionic threshold functions, respectively. One ob- 
serves that for k 3> T both threshold functions essentially behave as for zero temperature. 

4 For the present choice of Rk the temperature dependence of the threshold functions is considerably 
smoother than in that of previous investigations fl32j . 
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Figure 7: The plot shows the temperature dependence of the bosonic (a) and the fermionic (b) threshold 
functions lf{w, T/k) and l[ F ^ 4 (w, T/k), respectively, for different values of the dimensionless mass term w. 
The solid line corresponds to w = whereas the dotted ones correspond to w = 0.1, w = 1 and w = 10 
with decreasing size of the dots. For T 3> k the bosonic threshold function becomes proportional to T/k 
whereas the fermionic one tends to zero. In this range the theory with properly rescaled variables behaves 
as a classical three-dimensional theory. 



38 



For growing T or decreasing k this changes as more and more Matsubara modes decouple 
until finally all massive modes are suppressed. The bosonic threshold function If shows for 
k ^ T the linear dependence on T/k derived in Eq.(10.4). In particular, for the bosonic 



excitations the threshold function for w < 1 can be approximated with reasonable accu- 
racy by Z^(u>;?7$) for T/k < 0.25 and by (4T/fc)/^(w; r/$) for T/k > 0.25. The fermionic 
threshold function tends to zero for k <C T since there is no massless fermionic zero 
mode, i.e. in this limit all fermionic contributions to the flow equations are suppressed. On 
the other hand, the fermions remain quantitatively relevant up to T/k ~ 0.6 because of 
the relatively long tail in fig. ^b. The transition from four to three-dimensional threshold 
functions leads to a smooth dimensional reduction as k is lowered from k ^> T to k <C T\ 
Whereas for k 3> T the model is most efficiently described in terms of standard four- 
dimensional fields $ a choice of rescaled three-dimensional variables $3 = $/ a/T becomes 
better adapted for k < T. Accordingly, for high temperatures one will use a potential 

u S (t,p 3 ) = ^u(t,p) ; p 3 = . (10.5) 

In this regime corresponds to the free energy of classical statistics and Tk>o is a 

classical coarse grained free energy. 

For our numerical calculations at non-vanishing temperature we exploit the discussed 
behavior of the threshold functions by using the zero temperature flow equations in the 
range k > 10T. For smaller values of k we approximate the infinite Matsubara sums 
(cf. Eq.( |10.3j )) by a finite series such that the numerical uncertainty at k = 10T is better 
than 10~ 4 . This approximation becomes exact in the limit k <C 10T. 



11 The quark meson model at T ^ 

So far we have considered the relevant fluctuations that contribute to the flow of in 
dependence on the scale k. In a thermal equilibrium situation also depends on the 
temperature T and one may ask for the relevance of thermal fluctuations at a given scale 
k. In particular, for not too high values of T the "initial condition" for the solution of 
the flow equations should essentially be independent of temperature. This will allow us to 
fix from phenomenological input at T = and to compute the temperature dependent 
quantities in the infrared (k — > 0). We note that the thermal fluctuations which contribute 
to the r.h.s. of the flow equation for the meson potential Eq.( |8.2| ) are effectively suppressed 
for T ^ k/4 as discussed in detail in the last section. Clearly for T <: /c$/3 temperature 
effects become important at the compositeness scale. We expect the linear quark meson 
model with a compositeness scale k$ ~ 600 MeV to be a valid description for two flavor 
QCD below a temperature of about 170 MeV. We note that there will be an effective 
temperature dependence of induced by the fluctuations of other degrees of freedom 
besides the quarks, the pions and the sigma which are taken into account here. We will 
comment on this issue in the last section. For realistic three flavor QCD the thermal kaon 
fluctuations will become important for T <; 170 MeV. 
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Figure 8: The expectation value 2cto is shown as a function of temperature T for three different values 
of the zero temperature pion mass. 



We compute the quantities of interest for temperatures T ^ 170 MeV by numerically 
solving the T-dependent version of the flow equations (|8.2j) , ( |3.5|) — ( ^.7|) by lowering 
k from kq> to zero. For this range of temperatures we use the initial values as given in 
the first line of tab. This corresponds to choosing the zero temperature pion mass and 
the pion decay constant (f n = 92.4 MeV for m n = 135 MeV) as phenomenological input. 
The only further input is the constituent quark mass M q which we vary in the range 
M q ~ 300 — 350 MeV. We observe only a minor dependence of our results on M q for the 
considered range of values. In particular, the value for the critical temperature T c of the 
model remains almost unaffected by this variation. 

We have plotted in fig. |8| the renormalized expectation value 2<7o of the scalar field as a 
function of temperature for three different values of the average light current quark mass 
rh. (We remind that 2a"o(T = 0) = /„-.) For rh = the order parameter o"o of chiral 
symmetry breaking continuously goes to zero for T — > T c = 100.7 MeV characterizing the 
phase transition to be of second order. The universal behavior of the model for small T — T c 
and small rh is discussed in more detail in the following section. We point out that the value 
of T c corresponds to = 80.8 MeV, i.e. the value of the pion decay constant for rh — 0, 
which is significantly lower than / w = 92.4 MeV obtained for the realistic value m p h y s- If 
we would fix the value of the pion decay constant to be 92.4 MeV also in the chiral limit 
(m = 0), the value for the critical temperature would raise to 115 MeV. The nature of the 
transition changes qualitatively for m/fl where the second order transition is replaced by a 
smooth crossover. The crossover for a realistic m p h ys or m 7r (T = 0) = 135 MeV takes place 
in a temperature range T ~ (120 — 150) MeV. The middle curve in fig. |8| corresponds to a 
value of rh which is only a tenth of the physical value, leading to a zero temperature pion 
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Figure 9: The plot shows the m a as a function of temperature T for three different values of the zero 
temperature pion mass. 



mass m n = 45 MeV. Here the crossover becomes considerably sharper but there remain 
substantial deviations from the chiral limit even for such small quark masses rh ~ 1 MeV. 
The temperature dependence of has already been mentioned (see fig. |6|) for the same 
three values of rh. As expected, the pions behave like true Goldstone bosons for rh = 0, 
i.e. their mass vanishes for T < T c . Interestingly, m n remains almost constant as a function 
of T for T < T c before it starts to increase monotonically. We therefore find for two flavors 
no indication for a substantial decrease of m v around the critical temperature. 

The dependence of the mass of the sigma resonance m a on the temperature is displayed 
in fig. [| for the above three values of rh. In the absence of explicit chiral symmetry break- 
ing, rh = 0, the sigma mass vanishes for T < T c . For T < T c this is a consequence of the 
presence of massless Goldstone bosons in the Higgs phase which drive the renormalized 
quartic coupling A to zero. In fact, A runs linearly with k for T ^ k/4 and only evolves log- 
arithmically for T ^ k/4. Once m^fl the pions acquire a mass even in the spontaneously 
broken phase and the evolution of A with k is effectively stopped at k <~ m n . Because of 
the temperature dependence of 00^=0 (cf. fig. H) this leads to a monotonically decreasing 
behavior of m a with T for T ^ T c . This changes into the expected monotonic growth 
once the system reaches the symmetric phase for T > T c . For low enough rh one may use 
the minimum of m cr (T) for an alternative definition of the (pseudo-) critical temperature 
denoted as Tj£- Tab. in the introduction shows our results for the pseudocritical tem- 
perature for different values of rh or, equivalently, m n (T = 0). For a zero temperature pion 
mass m n = 135 MeV we find Tj£' = 128 MeV. At larger pion masses of about 230 MeV 
we observe no longer a characteristic minimum for m a apart from a very broad, slight dip 
at T ~ 90 MeV. A comparison of our results with lattice data is given below in the next 
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Figure 10: The plot shows the derivative of the meson potential U(T) with respect to the renormalizcd 
field 4>r = for different values of T. The first curve on the left corresponds to T — 175 MeV. 

The successive curves to the right differ in temperature by AT = 10 MeV down to T = 5 MeV. 



section. 

Our results for the chiral condensate (^>ip\ as a function of temperature for different val- 
ues of the average current quark mass are presented in fig. |5|. We will compare l^if) (T, rh) 
with its universal scaling form for the 0(4) Heisenberg model in the following section. 

Our ability to compute the complete temperature dependent effective meson potential 



U is demonstrated in fig. |10| where we display the derivative of the potential with respect 
to the renormalized field 4>r = , for different values of T. The curves cover a 

temperature range T = (5—175) MeV. The first one to the left corresponds to T = 175 MeV 
and neighboring curves differ in temperature by AT = 10 MeV. One observes only a weak 
dependence of dU{T)/d4>R on the temperature for T ^ 60 MeV. Evaluated for (p R = a 
this function connects the renormalized field expectation value with m n (T), the source j 
and the mesonic wave function renormalization Z<p(T) according to 

= *o) = = 4a (TK( T) . (1L1) 

We close this section with a short assessment of the validity of our effective quark 
meson model as an effective description of two flavor QCD at non-vanishing temperature. 
The identification of qualitatively different scale intervals which appear in the context of 
chiral symmetry breaking, as presented in the preceding sections for the zero temperature 
case, can be generalized to T ^ 0: For scales below kq> there exists a hybrid description in 
terms of quarks and mesons. For k x $B < k ^ 600 MeV chiral symmetry remains unbroken 
where the symmetry breaking scale k x ss{T) decreases with increasing temperature. Also 
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the constituent quark mass decreases with T. The running Yukawa coupling depends only 
mildly on temperature for T ^ 120 MeV. (Only near the critical temperature and for 
m = the running is extended because of massless pion fluctuations.) On the other hand, 
for k ^ AT the effective three-dimensional gauge coupling increases faster than at T = 
leading [|T^] to an increase of Aqcd(T) with T. As k gets closer to the scale Aqcd(T) it is 
no longer justified to neglect in the quark sector confinement effects which go beyond the 
dynamics of our present quark meson model. Here it is important to note that the quarks 
remain quantitatively relevant for the evolution of the meson degrees of freedom only for 
scales k ^ T/0.6 (cf. fig. 0). In the limit k <C T/0.6 all fermionic Matsubara modes decouple 
from the evolution of the meson potential according to the temperature dependent version 
of Eq. ( |3.2|) . Possible sizeable confinement corrections to the meson physics may occur if 
Aqcd(^) becomes larger than the maximum of M q (T) and T/0.6. This is particularly 
dangerous for small m in a temperature interval around T c . Nevertheless, the situation 
is not dramatically different from the zero temperature case since only a relatively small 
range of k is concerned. We do not expect that the neglected QCD non-localities lead to 
qualitative changes. Quantitative modifications, especially for small m and \T — T c \ remain 
possible. This would only effect the non-universal amplitudes which will be discussed in 
the next section. The size of these corrections depends on the strength of (non-local) 
deviations of the quark propagator and the Yukawa coupling from the values computed in 
the quark meson model. 



12 Critical behavior near the chiral phase transition 

In this section we study the linear quark meson model in the vicinity of the critical tem- 
perature T c close to the chiral limit rh = 0. In this region we find that the sigma mass m~ l 
is much larger than the inverse temperature T -1 , and one observes an effectively three- 
dimensional behavior of the high temperature quantum field theory. We also note that the 
fermions are no longer present in the dimensionally reduced system as has been discussed 
above. We therefore have to deal with a purely bosonic 0(4)-symmetric linear sigma 
model. At the phase transition the correlation length becomes infinite and the effective 
three-dimensional theory is dominated by classical statistical fluctuations. In particular, 
the critical exponents which describe the singular behavior of various quantities near the 
second order phase transition are those of the corresponding classical system. 

Many properties of this system are universal, i.e. they only depend on its symmetry 
(0(4)), the dimensionality of space (three) and its degrees of freedom (four real scalar 
components). Universality means that the long-range properties of the system do not 
depend on the details of the specific model like its short distance interactions. Neverthe- 
less, important properties as the value of the critical temperature are non-universal. We 
emphasize that although we have to deal with an effectively three-dimensional bosonic 
theory, the non-universal properties of the system crucially depend on the details of the 
four-dimensional theory and, in particular, on the fermions. 

Our aim is a computation of the critical equation of state which relates for arbitrary 
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T near T c the derivative of the free energy or effective potential U to the average current 
quark mass rh. The equation of state then permits to study the temperature and quark 
mass dependence of properties of the chiral phase transition. 

At the critical temperature and in the chiral limit there is no scale present in the theory. 
In the vicinity of T c and for small enough rh one therefore expects a scaling behavior of the 
effective average potential and accordingly a universal scaling form of the equation of 



state [J28[] . There are only two independent scales close to the transition point which can be 
related to the deviation from the critical temperature, T — T c , and to the explicit symmetry 
breaking through the quark mass rh. As a consequence, the properly rescaled potential 
can only depend on one scaling variable. A possible choice for the parameterization of the 
rescaled "unrenormalized" potential is the use of the Widom scaling variable 



(T - T c ) IT C 

Here f3 is the critical exponent of the order parameter a in the chiral limit rh = (see 
Eq.( |12.5|) ). With U'(p = 2<7q) = j/(2<t ) the Widom scaling form of the equation of state 



reads [^BJ 



-i- ; , , , f(x) (12.2) 

where the exponent 5 is related to the behavior of the order parameter according to 
Eq.( 12.7 ). The equation of state Eq. ([L2.2 ) is written for convenience directly in terms 



of four-dimensional quantities. They are related to the corresponding effective variables 
of the three-dimensional theory by appropriate powers of T c . The source j is determined 
by the average current quark mass rh as j = 2m\ i> rh. The mass term at the composite- 
ness scale, fnl^, also relates the chiral condensate to the order parameter according to 

(ipip^ = — 2m| <i) (o r o — rh). The critical temperature of the linear quark meson model was 
found above to be T c = 100.7 MeV. 

The scaling function / is universal up to the model specific normalization of x and 
itself. Accordingly, all models in the same universality class can be related by a rescaling 
of Zfo and T — T c . The non-universal normalizations for the quark meson model discussed 
here are defined according to 

f(0) = D , f{-B-W) = . (12.3) 

We find D = 1.82 ■ 10 -4 , B = 7.41 and our result for (3 is given in tab. |||. Apart from 
the immediate vicinity of the zero of f(x) we find the following two parameter fit for the 
scaling function |2!| , 

/a (a) = 1.816 -10" 4 (1 + 136.1a;) 2 (1 + 160.9 9x) A 

(1 + 160.9 (0.9446 ^-i/( 7 -2-A) ^ 7 - 2 -A ( l 2 ' l) 

to reproduce the numerical results for / and df jdx at the 1—2% level with 9 = 0.625 (0.656), 
A = —0.490 (—0.550) for x > (x < 0) and 7 as given in tab. |3|. The universal properties 
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of the scaling function can be compared with results obtained by other methods for the 
three-dimensional 0(4) Heisenberg model. In fig. |ll] we display our results along with those 



obtained from lattice Monte Carlo simulation |p7 |, second order epsilon expansion |)8| and 



mean field theory. We observe a good agreement of average action, lattice and epsilon 
expansion results within a few per cent for T < T c . Above T c the average action and 
the lattice curve go quite close to each other with a substantial deviation from the epsilon 
expansion and mean field scaling function. (We note that the question of a better agreement 
of the curves for T < T c or T > T c depends on the chosen non-universal normalization 
conditions for x and / (cf. Eq. (|12.3| )).) 

Before we use the scaling function f(x) to discuss the general temperature and quark 
mass dependent case, we consider the limits T = T c and rh = 0, respectively. In these 
limits the behavior of the various quantities is determined solely by critical amplitudes and 
exponents. In the spontaneously broken phase (T < T c ) and in the chiral limit we observe 
that the renormalized and unrenormalized order parameters scale according to 

2 °°(T) , ocV/2 (T C -TY' 2 



( 2E Y , r 

; ' (i2.5) 

2<x (T) fTc -T\ p 



T c V T c 

respectively, with E = 0.814 and the value of B given above. In the symmetric phase the 
renormalized mass m = m n = m a and the unrenormalized mass m = Z^ 2 m behave as 

mm - {e r ( T -^X , 

/, (12.6) 
m(T) ^ +v i/2 ( T-T c y' 2 ^ 



T r ^ ) V T, 



T r V / V T, 



- C 



where £ + = 0.270, C + = 2.79. For T = T c and non- vanishing current quark mass we have 

(12.7) 



2a < ^ ^ 1/5 



T, 



c 




with the value of D given above. 

Though the five amplitudes E, B, £ + , C + and D are not universal there are ratios of 
amplitudes which are invariant under a rescaling of Zfo and T — T c . Our results for the 
universal amplitude ratios are 



R x 


= C + DB 




= {tf /u 


tE 


= 0.220. 



V+VB = 0.852 , (12.8) 
Those for the critical exponents are given in tab. [3]. Here the value of 77 is obtained 
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Figure 11: The figure shows a comparison of our results, denoted by "average action", with results 
of other methods for the scaling function of the three-dimensional 0(4) Heisenberg model. We have 
labeled the axes for convenience in terms of the expectation value (To and the source j of the corresponding 
four-di mensi onal theory. The constants B and D specify the non-universal amplitudes of the model 
(cf. Eq.( 12.3 ). The curve labeled by "MC" represents a fit to lattice Monte Carlo data. The second order 
epsilon expansion and mean field results are denoted by "e" and "mf" , respectively. Apart from our 
results the curves are taken from |67|1. 
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V 


7 


5 


P 


V 


average action 


0.787 


1.548 


4.80 


0.407 


0.0344 


3d PT 


0.73(2) 


1.44(4) 


4.82(5) 


0.38(1) 


0.03(1) 


MC 


0.7479(90) 


1.477(18) 


4.851(22) 


0.3836(46) 


0.0254(38) 



Table 3: The table shows the critical exponents corresponding to the three-dimensional 0(4)-Heisenberg 
model. Our results are denoted by "average action" whereas "3d PT" labels the exponents obtained from 
a 3d perturbative expansion using Pade techniques J7(], fFl| . The bottom line contains lattice Monte Carlo 
results. J6?J 



from the temperature dependent version of Eq. (|8.7j ) at the critical temperature |52|| . For 



comparison tab. also gives the results |7(], |7T|] of a 3d perturbative computation as well 
as lattice Monte Carlo results []69| which have been used for the lattice form of the scaling 
function in fig. [11]. f] There are only two independent amplitudes and critical exponents, 
respectively. They are related by the usual scaling relations of the three-dimensional scalar 
0(iV)-model |7T[ which we have explicitly verified by the independent calculation of our 
exponents. 

We turn to the discussion of the scaling behavior of the chiral condensate (vpip^ for the 
general case of a temperature and quark mass dependence. In fig. || we have displayed our 
results for the scaling equation of state in terms of the chiral condensate 




as a function of T/T c = 1 + x(j/T^f(x)Y^ s for different quark masses or, equivalently, 
different values of j. The curves shown in fig. correspond to quark masses m = 0, 
m = m p hys/10, m = m p h ys and m = 3.5m p h ys or, equivalently, to zero temperature pion 
masses m w = 0, m w = 45 MeV, m w = 135 MeV and = 230 MeV, respectively (cf. fig. |). 
One observes that the second order phase transition with a vanishing order parameter at 
T c for m = is turned into a smooth crossover in the presence of non-zero quark masses. 
The scaling form Eq. (|12.9|) for the chiral condensate is exact only in the limit T — > T ( 



CI 



j — ► 0. It is interesting to find the range of temperatures and quark masses for which (i/iip 
approximately shows the scaling behavior Eq.( |12.9|) . This can be infered from a comparison 



(see fig. |5|) with our full non-universal solution for the T and j dependence of ypip 
For m n = one observes approximate scaling behavior for temperatures T ^ 90 MeV. 
This situation persists up to a pion mass of m n = 45 MeV. Even for the realistic case, 
m v = 135 MeV, and to a somewhat lesser extent for = 230 MeV the scaling curve 
reasonably reflects the physical behavior for T ^ T c . For temperatures below T c , however, 



5 See also ref ]7^] and references therein for a recent calculation of critical exponents using similar 
methods as in this work. For high precision estimates of the critical exponents see also refs |73, Q . 



47 



the zero temperature mass scales become important and the scaling arguments leading to 
universality break down. 

The above comparison may help to shed some light on the use of universality arguments 
away from the critical temperature and the chiral limit. One observes that for temperatures 
above T c the scaling assumption leads to quantitatively reasonable results even for a pion 
mass almost twice as large as the physical value. This in turn has been used for two 
flavor lattice QCD as theoretical input to guide extrapolation of results to light current 
quark masses. From simulations based on a range of pion masses 0.3 ^ m^/nip ^ 0.7 
and temperatures < T ^ 250 MeV a "pseudocritical temperature" of approximately 
140 MeV with a weak quark mass dependence is reported Jn| . Here the "pseudocritical 
temperature" T pc is defined as the inflection point of (vpip^ as a function of temperature. 
The values of the lattice action parameters used in |75|] with N t = 6 were am = 0.0125, 
Q/g 2 = 5.415 and am = 0.025, Q/g 2 = 5.445. For comparison with lattice data we have 
displayed in fig. ^| the temperature dependence of the chiral condensate for a pion mass 
= 230 MeV. From the free energy of the linear quark meson model we obtain in 
this case a pseudocritical temperature of about 150 MeV in reasonable agreement with 
the results of ref ||75|| . In contrast, for the critical temperature in the chiral limit we 



obtain T c = 100.7 MeV. This value is considerably smaller than the lattice results of about 
(140—150) MeV obtained by extrapolating to zero quark mass in ref |f75| . We point out that 
for pion masses as large as 230 MeV the condensate (ipip^ (T) is almost linear around the 
inflection point for quite a large range of temperature. This makes a precise determination 
of T c somewhat difficult. Furthermore, fig. [5] shows that the scaling form of \M)ip) (T) 
underestimates the slope of the physical curve. Used as a fit with T c as a parameter this 
can lead to an overestimate of the pseudocritical temperature in the chiral limit. We also 
mention here the results of ref |76|| . There two values of the pseudocritical temperature, 
T pc = 150(9) MeV and T pc = 140(8), corresponding to am = 0.0125, Q/g 2 = 5.54(2) and 
am = 0.00625, 6/g 2 = 5.49(2), respectively, (both for N t = 8) were computed. These 
values show a somewhat stronger quark mass dependence of T pc and were used for a linear 
extrapolation to the chiral limit yielding T c = 128(9) MeV. 

The linear quark meson model exhibits a second order phase transition for two quark 
flavors in the chiral limit. As a consequence the model predicts a scaling behavior near 
the critical temperature and the chiral limit which can, in principle, be tested in lattice 
simulations. For the quark masses used in the present lattice studies the order and univer- 
sality class of the transition in two flavor QCD remain a partially open question. Though 
there are results from the lattice giving support for critical scaling ]77|, |78[ there are also 
recent simulations with two flavors that reveal significant finite size effects and problems 



with 0(4) scaling [79, 8C 



13 Additional degrees of freedom 

So far we have investigated the chiral phase transition of QCD as described by the linear 
0(4)-model containing the three pions and the sigma resonance as well as the up and 
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down quarks as degrees of freedom. Of course, it is clear that the spectrum of QCD is 
much richer than the states incorporated in our model. It is therefore important to ask to 
what extent the neglected degrees of freedom like the strange quark, strange (pseudo) scalar 
mesons, (axial) vector mesons, baryons, etc., might be important for the chiral dynamics 
of QCD. Before doing so it is perhaps instructive to first look into the opposite direction 
and investigate the difference between the linear quark meson model described here and 



chiral perturbation theory based on the non-linear sigma model |36[] . In some sense, chiral 
perturbation theory is the minimal model of chiral symmetry breaking containing only the 
Goldstone degrees of freedom. By construction it is therefore only valid in the sponta- 
neously broken phase and can not be expected to yield realistic results for temperatures 
close to T c or for the symmetric phase. However, for small temperatures (and momentum 
scales) the non-linear model is expected to describe the low-energy and low-temperature 
limit of QCD reliably as an expansion in powers of the light quark masses. For vanish- 



ing temperature it has been demonstrated recently |43], [0J that the results of chiral 
perturbation theory can be reproduced within the linear meson model once certain higher 
dimensional operators in its effective action are taken into account for the three flavor 
case. Moreover, some of the parameters of chiral perturbation theory (L 4 , . . . , L 8 ) can be 
expressed and therefore also numerically computed in terms of those of the linear model. 
For non-vanishing temperature one expects agreement only for low T whereas deviations 
from chiral perturbation theory should become large close to T c . Yet, even for T <^ T c 
small quantitative deviations should exist because of the contributions of (constituent) 
quark and sigma meson fluctuations in the linear model which are not taken into account 
in chiral perturbation theory. 

From [[81] we infer the three-loop result for the temperature dependence of the chiral 



condensate in the chiral limit for N light flavors 
(^)CZVr = (^) xPr (0) 
+ iV(iV 2 -l) 




:i3.i} 



The scale Ti can be determined from the .D-wave isospin zero tttt scattering length and is 
given by Ti = (470 ±100) MeV. The constant F Q is (in the chiral limit) identical to the pion 
decay constant F = = 80.8 MeV. In fig. [12] we have plotted the chiral condensate 



as a function of T/F for both, chiral perturbation theory according to Eq.( 13.1 ) and 
for the linear quark meson model. As expected the agreement for small T is very good. 
Nevertheless, the anticipated small numerical deviations present even for T < T c due to 
quark and sigma meson loop contributions are manifest. For larger values of T, say for 
T 0.8f^ the deviations become significant because of the intrinsic inability of chiral 
perturbation theory to correctly reproduce the critical behavior of the system near its 
second order phase transition. 

Within the language of chiral perturbation theory the neglected effects of thermal quark 
fluctuations may be described by an effective temperature dependence of the parameter 
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T/fi 0) 

Figure 12: The plot displays the chiral condensate {ipip) as a function of T/fi°\ The solid line corre- 
sponds to our results for vanishing average current quark mass rh = whereas the dashed line shows the 
corresponding three-loop chiral perturbation theory result for I\ = 470 MeV. 

F (T). We notice that the temperature at which these corrections become important equals 
approximately one third of the constituent quark mass M q (T) or the sigma mass m a (T), 
respectively, in perfect agreement with fig. 0. As suggested by this figure the onset of the 
effects from thermal fluctuations of heavy particles with a T-dependent mass m^(T) is 
rather sudden for T ^ m#(T)/3. These considerations also apply to our two flavor quark 
meson model. Within full QCD we expect temperature dependent initial values at k$. 

The dominant contribution to the temperature dependence of the initial values presum- 
ably arises from the influence of the mesons containing strange quarks as well as the strange 
quark itself. Here the quantity m? seems to be the most important one. (The temperature 
dependence of higher couplings like A(T) is not very relevant if the IR attractive behavior 
remains valid, i.e. if Z<$,^$ remains small for the range of temperatures considered. We 
neglect a possible T-dependence of the current quark mass m.) In particular, for three 
flavors the potential contains a term 

- ^I7 fc$ (det $ + det $ f ) = -V ktb ip s $ uu § dd + ... (13.2) 

which reflects the axial £7a(1) anomaly. It yields a contribution to the effective mass term 
proportional to the expectation value ($ ss ) = <f s , i.e. 

= ~\ V k^s ■ (13.3) 

Both, V kii and ip s , depend on T. We expect these corrections to become relevant only for 
temperatures exceeding m^(T)/3 or M s (T)/3. We note that the temperature dependent 
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kaon and strange quark masses, mxiT) and M S (T), respectively, may be somewhat differ- 
ent from their zero temperature values but we do not expect them to be much smaller. A 
typical value for these scales is around 500 MeV. Correspondingly, the thermal fluctuations 
neglected in our model should become important for T ^ 170 MeV. It is even conceivable 
that a discontinuity appears in tp s {T) for sufficiently high T (say T ~ 170 MeV). This 
would be reflected by a discontinuity in the initial values of the 0(4)-model leading to a 
first order transition within this model. Obviously, these questions should be addressed 
in the framework of the three flavor SUl{3) x SUr(3) quark meson model. Work in this 
direction is in progress. 

We note that the temperature dependence of U{T)(p s (T) is closely related to the question 
of an effective high temperature restoration of the axial C/a(1) symmetry [|58, The rj 
mass term is directly proportional to this combination pE2| , m^,(T)—ml(T) ~ ^u(T)(p s (T). 
Approximate Z7a(1) restoration would occur if (f s {T) or v(T) would decrease sizeable for 
large T. For realistic QCD this question should be addressed by a three flavor study. 
Within two flavor QCD the combination V^ s is replaced by an effective anomalous mass 
term . The temperature dependence of (T) could be studied by introducing quarks 



and the axial anomaly in the two flavor matrix model of ref [54|. We add that this question 



has also been studied within full two flavor QCD in lattice simulations [79, p3, B3j. So far 



there does not seem to be much evidence for a restoration of the Ua{X) symmetry near T c 
but no final conclusion can be drawn yet. 

To summarize, we have found that the effective two flavor quark meson model presum- 
ably gives a good description of the temperature effects in two flavor QCD for a temperature 
range T ^ 170 MeV. Its reliability should be best for low temperature where our results 
agree with chiral perturbation theory. However, the range of validity is considerably ex- 
tended as compared to chiral perturbation theory and includes, in particular, the critical 
temperature of the second order phase transition in the chiral limit. We have explicitly 
connected the universal critical behavior for small \T — T c \ and small current quark masses 
with the renormalized couplings at T = and realistic quark masses. The main quanti- 
tative uncertainties from neglected fluctuations presumably concern the values of and 
T c which, in turn, influence the non-universal amplitudes B and D in the critical region. 
We believe that our overall picture is rather solid. Where applicable our results compare 
well with numerical simulations of full two flavor QCD. 



14 Conclusions 

Our conclusions may be summarized in the following points: 

1. The connection between short distance perturbative QCD and long distance meson 
physics by analytical methods seems to emerge step by step. The essential ingredients 
are nonperturbative flow equations as approximations of exact renormalization group 
equations and a formalism which allows a change of variables by introducing meson- 
like composite fields. 
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2. The relevance of meson-like ipip composite objects is established in this framework. 
A typical scale where the mesonic bound states appear is the compositeness scale 
kq, ~ (600 — 700) MeV. The occurrence of chiral symmetry breaking depends on the 
effective meson mass at the compositeness scale. For a certain range of values of 
the ratio fnl^/k^, spontaneous chiral symmetry breaking is induced by quark fluctu- 
ations. A definite analytical establishment of spontaneous chiral symmetry breaking 
from "first principles" (i.e., short distance QCD) still awaits a reliable calculation of 
this ratio. 

3. Phenomenologically, the ratio ml^/k^ may be determined from the value of the 
constituent quark mass in units of the pion decay constant f n . For this value the 
ratios f n /k^ and (^pip^j /k% can be computed. Together with an earlier estimate of /c$ 

this yields rather encouraging results: f w ~ 100 MeV and (ipip^ — —(190 MeV) 3 for 
two flavor QCD. It will be very interesting to generalize these results to the realistic 
three flavor case. 

4. The formalism presented here naturally leads to an effective linear quark meson model 
for the description of mesons below the compositeness scale k$. For this model the 
standard non-linear sigma model of chiral perturbation theory emerges as a low 
energy approximation due to the large sigma mass. 

5. The old puzzle about the precise connection between the current and the constituent 
quark mass reveals new interesting aspects in this formalism. In the context of the 
effective average action one may define these masses by the quark propagator at zero 
or very small momentum, either in the quark gluon picture (current quark mass) or 
the effective quark meson model (constituent quark mass). There is no conceptual 
difference between the two situations. As a function of k the running quark mass 
smoothly interpolates between the standard current quark mass m q for high k and 
the standard constituent quark mass M q for low k. (This holds at least as long as 
the minimum of the effective scalar potential is unique.) A rapid quantitative change 
occurs for k ~ (400 — 500) MeV because of the onset of chiral symmetry breaking. 
For k = one expects this behavior to carry over to the momentum dependence of 
the quark propagator: For small q 2 the inverse quark propagator is dominated by 
M q . In contrast, for high q 2 the constant term in the inverse propagator is reduced 
to m q since q 2 replaces k 2 as an effective infrared cutoff. One expects a smooth 
interpolation between the two limits and it would be interesting to know the form 
of the propagator for momenta in the transition region. Furthermore, the symmetry 
breaking source term which determines the pion mass in the linear or non-linear 
meson model can be related to the current quark mass at the compositeness scale. 
This constitutes a bridge between the low energy meson properties and the running 
quark mass at a scale which is not too far from the validity of perturbation theory. 

6. A particular version of the Nambu-Jona-Lasinio model appears in our formalism as 
a limiting case (infinitely strong renormalized Yukawa coupling at the compositeness 
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scale k$). Here the ultraviolet cutoff which is implicit in the NJL model is dictated 
by the momentum dependence of the infrared cutoff function in the effective 
average action. The characteristic cutoff scale is k$. In this context our results can 
be interpreted as an approximative solution of the NJL model. Our method includes 
many contributions beyond the leading order contribution of the 1/N C expansion. 

7. Based on a satisfactory understanding of the meson properties in the vacuum we 
have described their behavior in a thermal equilibrium situation. Our method should 
remain valid for temperature below ~ 170 MeV. This extends well beyond the validity 
of chiral perturbation theory. In particular, in the chiral limit of vanishing quark 
masses we can describe the universal critical behavior near a second order phase 
transition of two flavor QCD. The universal behavior is quantitatively connected to 
observed quantities at zero temperature and realistic quark masses. 
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